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1 . I ntroduction 

A variety  of  methods  is  possible  to  attack  hydrodynamic 
problems.  Particularly,  two  different  methods  are  most 
widely  used  to  obtain  useful  information  for  the  flow 
characteristics  around  a hydrofoil  or  propeller  blades. 

One  is  the  series  expansion  methods  of  a small  parameter. 
This  type  of  method  is  known  as  a singular  perturbation 
problem  or  "the  method  of  matched  asymptotic  expansions." 

The  second  method  is  the  purely  numerical  one.  Undoubtedly, 
i his  method  has  been  becoming  of  increasing  importance  with 
he  development  of  modern  computing  devices.  Therefore, 
nr  thods  employing  numerical  lifting  surface  theory  are  now 
generally  considered  as  essential  for  the  satisfactory 
d>-.,ign  and  performance  prediction  of  marine  propellers. 

Since  1960,  a systematic  research  on  these  subjects  has 
b'  en  carried  out  at  MIT  under  the  direction  of  Professor 
kerwin  and,  as  a result,  a series  of  computer  programs  has 
been  developed  such  as  MIT  PBD  , PINV  and  PUF  -series. 

T i newest  program,  PUF-2,  which  can  be  applied  to  both 
toady  and  unsteady  cases  to  determine  the  hydrodynamic 

* Designation  of  MIT  "propeller  blade  section  design", 
"propeller  inverse"  and  "propeller  unsteady  force" 
programs,  respectively. 
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loading  on  a propeller  of  prescribed  geometry,  is  believed 
to  be  one  of  the  best  presently  available  numerical  schemes 
and  to  offer  the  greatest  potential  for  the  determination 
of  detailed  hydrodynamic  loading  around  blades.  In  order 
to  verify  the  overall  accuracy  and  achieve  the  possible 
further  refinement  of  the  numerical  scheme,  however,  it  is 
necessary  to  conduct  a series  of  appropriate  experiments. 

Until  very  recently,  the  experimental  confirmation  of 
numerical  lifting  surface  theory  has  been  limited  to  a 
comparison  of  overall  propeller  blade  thrust  and  torque 
loading.  For  many  years,  experimental  fluid  mechanics  has 
made  use  of  standard  mechanical  measuring  probes  to  obtain 
information  on  fluid  velocity.  Total-pressure  probes,  in 
conjunction  with  static-pressure  probes,  have  provided  the 
principal  means  of  measuring  mean  velocity  and  hot-wire  or 
hot-film  anemometers  the  principal  means  of  measuring  instan- 
taneous velocity  and,  hence,  mean  velocity,  rms  velocities 
and  velocity  correlations.  Mechanical  probes  will  undoubt- 
edly continue  to  be  important  to  experimental  fluid  mechanics, 
but  it  is  obvious  that  an  experiment  for  the  detailed  load 
distribution  around  a propeller  blade  with  such  conventional 
measuring  devices  is  extremely  difficult  and  costly,  if  not 
impossible.  The  limitations  on  the  conventional  mechanical 
probes  make  new  measuring  techniques  worthy  of  careful  con- 
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sideration.  Thus,  "Laser-Doppler  Anomometry  Technique" 

has  been  introduced  as  an  experimental  method  expecting 

to  be  able  to  overcome  all  of  the  limitations  that  the 

mechanical  probes  have.  At  the  stage  of  paper  survey, 

the  author  had  encountered  several  different  terminologies 

used  from  paper  to  paper  regarding  this  technique  such 

as  laser-Doppler  anemometry,  laser  anemometry,  optical 

anemometry,  laser-Doppler  velocimetry  £tc.  However,  the 

term  "Laser-Doppler  Anemometry",  or  simply  "LDA"  will  be 

used  exclusively  in  this  paper.  The  strong  interest  in 

LDA  is  evidenced  by  over  thousand  papers  that  have  been 

published  on  the  subject  since  the  first  one  by  Yeh  and 
• • __  * 

Cummings  in  1964  (_  1J  , who  measured  in  a fully-developed 
laminar  pipe  flow  of  water.  This  proliferation  of  papers 
is  due  both  to  the  potential  of  the  technique  as  well  as 
some  of  the  challenging  theoretical  considerations  asso- 
ciated with  it.  In  this  paper,  a brief  discussion  on  the 
basic  concept  of  the  LDA  technique,  necessary  components 
and  their  functions  will  first  be  introduced  along  with 
some  historical  review.  Then,  the  results  of  measurements 
and  comparison  with  the  numerical  output  will  be  discussed 
it  xt . 

* The  numbers  in  [ J denote  references  listed  at  the  end 
of  this  paper. 
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IT.  Experimental  Method 

11 . 1 . General 

As  mentioned  in  Chapter  I , LDA  technique  has  been 
exclusively  used  as  an  experimental  method.  In  section 
II. 2 » the  basic  concept  and  some  considerations  for  appli- 
cation will  be  discussed.  At  the  time  when  the  author 
first  began  to  employ  this  technique  in  the  flow  measure- 
ment situations ( late  1976),  he  was  unable  to  find  any  ref- 
erence systematically  describing  general  theory  and  appli- 
cation of  LDA  technique.  Therefore,  the  author  has  pre- 
pared this  material  for  the  purpose  of  providing  some 
basic  knowledge  about  this  technique  for  the  possible 
future  user.  In  fact,  this  chapter  is  not  only  the  sum- 
mary of  the  survey  of  presently  available  papers  in  this 
field,  but  also  the  result  from  the  author's  own  experience. 
Some  understanding  of  the  working  theory  would  be  helpful 
in  properly  applying  the  device  to  a particular  flow  problem. 

11. 2.  Laser-Doppler  Anemometry  Technique 

Laser  anemometry  is  a technique  which  utilizes  scattered 

4 

light  from  particulates  in  a fluid  to  measure  the  velocity  of 
that  fluid.  When  light  is  scattered  from  a moving  object,  a 
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stationary  observer  will  see  a change  in  the  frequency  of 
the  scattered  light  (Doppler  shift)  proportional  to  the 
velocity  of  the  object.  In  the  Laser-Dopplor  Anemometer, 
this  Doppler  shift  is  used  to  measure  the  velocity  of  par- 
ticles in  a fluid.  From  the  particle  velocity,  the  fluid 
velocity  is  inferred.  This  basic  concept  is  shown  schemat- 
ically in  Figure  1.  Theoretically,  any  kind  of  light  can 
be  used  as  an  illuminating  light  source.  However,  white 
light  such  as  room  light  or  day  light  has  too  large  spectral 
width  to  be  of  any  use  in  such  measurements.  In  practice,  a 
laser  beam  is  almost  exclusively  used  as  a light  source  due 
to  the  following  properties  of  the  laser  beam: 

1.  monochromatic 

2 . coherent 

3.  high  intensity 

4.  linearly  polarized 

5.  easy  focusing 

Fn  terms  of  its  potential,  the  advantage  which  make  the 
laser  anemometer  very  attractive  for  many  fluid  measuring 
appli cations  are: 


1.  No  probe  in  the  fluid. 

This  is  the  most  obvious  advantage.  Only  the  laser 
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liqht  itself  needs  to  enter  the  flow  field,  with  a window 
the  only  access.  This  can  be  vitally  important  in  flows 
where  a probe  would  cause  a disturbance  or  the  fluid  is 
hot,  corrosive,  or  otherwise  unsuitable  for  an  immersion 
of  a probe. 

2.  No  flow  calibration  required. 

In  theory,  the  laser  anemometer  is  linear.  It  will 
be  shown  later  that  the  laser  light  frequency,  fluid  re- 
fractive index  at  measuring  location,  and  optical  geometry 
uniquely  determine  the  relation( linear)  between  the  parti- 
cle velocity  and  the  frequency  output  of  the  photodetector. 
This  usually  eliminates  the  need  for  actual  velocity  cali- 

■k 

brations . 

3.  Precise  measurement  of  velocity  components. 

A single  component  LDA  system  measures  a single  com- 
ponent of  the  particle  velocity  and  is  independent  of  other 
velocity  components.  Adding  frequency  shifting  permits 
measurements  in  flows  when  the  velocity  reverses.  With 
a two  component  system,  the  individual  components  are 
measured  completely  independently. 


* The  overall  accuracy  is  affected  by  many  factors.  One  of 
them  is  the  precision  of  the  signal  processor ( frequency  to 
voltage  converter) . This  can  be  easily  checked  with  stan- 
dard electronic  test  instruments. 
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4 L 5 . Small  measuring  volume  and  high  frequency  response. 

These  are  essential  requirements  for  applications  i 

where  detailed  flow  measurements  are  required.  A rela- 
tively small  measuring  volume  and  inherently  fast  response 
gives  it  the  ability  to  follow  rapidly  changing  velocities 
in  the  fluid. 

The  above  advantages  must  be  quantif ied-particularly  the  ‘ 

last  two.  These  items  will  be  discussed  in  section  II. 3. 

A.(l).  These  characteristics  add  greatly  to  the  laser  ' 

anemometer's  value  as  a measurement  tool  and  have  made  it 
applicable  both  to  very  low  velocity  flows  as  well  as 
supersonic  flows. 

I 1.3.  Considerations  for  Selection  of  a System 

A complete  Laser  Doppler  Anemometer  includes,  at  the 
minimum,  the  following  components* 

1 . laser  source 

2.  transmitting  (focusing)  optics 

3.  receiving  (light  collecting)  optics 

4.  photodetector  to  convert  light  signal  to  eletrical 
signal 


5.  signal  processor  to  convert  frequency  to  a voltage 
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Other  components  are  added  to  accomplish  two  component 
measurements,  distinguish  flow  reversals,  provide  seeding, 
etc.  Several  parameters  must  be  considered  for  proper 
selection  and  application  of  a laser  anemometer  to  a flow 
situation.  Among  these  ares 

1.  the  optics  arrangement  or  "mode" 

2.  whether  seeding  the  flow  with  particles  is  reguired 
and,  if  it  is,  how  to  do  it 

3.  the  size  and  type  of  laser 

4.  the  type  of  signal  processor 

5.  the  type  of  data  analysis  eguipment  required 
Some  of  the  variables  that  affect  the  selection  ares 

1.  whether  only  one  or  more  than  one  velocity  component 
is  to  be  measured  at  one  time 

2.  whether  the  fluid  is  liquid  or  gas 

3.  the  expected  fluid  velocity  range 

4 c the  amount  and  type  of  natural  particulates  in  the 
fluid  and  the  allowable  concentration  of  particulates 

5,  the  accessibility  of  the  measuring  point 

6.  the  type  of  data  required,  i.e.,  mean  flow,  true  rms, 
spectrum,  correlations,  etc. 


Other  items,  not  included  above,  may  well  be  dominant  in  a 
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given  measurement  situation.  The  decision  on  the  above  items 
tend  to  be  interrelated  and  some  considerations  on  the  each 
of  items  will  be  discussed  in  the  following. 

A.  Selection  of  Optical  Arrangements  and  Lenses 

At  present,  the  following  three  basic  "modes"  are  used. 

1 . Reference  beam 

(local  oscillator  heterodyne  arrangement) 

2.  Dual  beam 

(fringe  system  or  differential  heterodyne  arrangement) 

3.  Dual  scatter 

(symmetric  heterodyne  arrangement) 

The  first  two  optical  arrangements  are  commonly  used,  but 
the  third  one  has  not  been  used  as  extensively  as  the  others. 
Most  early  measurements  were  made  with  the  "reference  beam" 
optical  mode.  In  this  mode,  the  laser  beam  is  split  into 
an  intense  scattering  beam  and  a weak  reference  beam.  The 
reference  beam  is  directed  on  to  a photocathode  where  it 
beats  with  light  scattered  from  the  strong  beam  by  particles 
moving  with  the  flow.  The  frequency  of  the  scattered  light 
has  been  altered  by  the  Doppler  effect  and  the  interference 
with  the  reference  beam  provides  a frequency  difference 
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which  is  directly  proportional  to  the  particle  velocity. 

This  arrangement  was  employed  in  the  pioneer  work  of  Yeh 
and  Cummings  and  has  subsequently  been  used  by  many  authors. 
Goldstein  and  Hagen  in  1967  Welch  and  Tomme  in  1967  [l~f 

and  Pike,  Jackson,  Bourke  and  Page  in  1968  are  considered 

as  early  persons  who  applied  it  to  turbulent  water  flows. 
This  system  is  often  useful  when  the  particle  concentration 
is  high  to  give  an  essentially  continuous  signal.  This  sim- 
plified the  signal  processing  problem  as  well  as  giving  an 
output  similar  to  a hot-wire.  However,  noise  due  to  multiple 
particles  in  the  measuring  volume  as  well  as  other  practical 
limitations  of  high  particle  concentration  led  to  development 
of  techniques  using  low  particle  concentrations.  The  "dual 
beam"  optical  mode  and  greatly  improved  signal  processors 
have  now  made  LDA  available  for  essentially  routine  measure- 
ments in  a wide  range  of  applications.  Therefore,  the  dual 
beam  mode  is  the  most  commonly  used  optical  arrangement 
nowadays.  One  of  the  earliest  application  of  this  mode 
may  be  the  example  by  Durst  and  Whitelaw  in  1971  to  measure 
mean  and  fluctuating  velocity  components  in  fully  developed 
channel  flow  of  water  [5,  6l.  Two  different  scattering 
methods  can  be  used  with  dual  beam  system.  Forward  scatter 
requires  the  least  laser  power  while  backscatter  requires 

only  one  window  and  is  functionally  more  convenient.  The 
obvious  advantage  of  dual  beam  backscatter  mode  is  the 
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rigidity  of  a system  where  the  entire  anemometer  is  bolted 
rigidly  together.  The  disadvantage  is  the  relatively  weak 
scattered  light  signal  when  compared  to  direct  forward 
scatter  (signal  strength  is  often  only  a percent  or  less 
of  the  forward  scatter  value.)  The  dual  scatter  mode  has 
found  its  primary  application  in  simultaneous  measurement 
of  two  flow  components  with  a single  laser.  Although  the 
small  acceptable  aperture  size  (due  to  coherence  require- 
ments) limits  the  signal  to  noise  ratio,  the  simplicity 
of  the  system  for  two  component  measurements  is  attractive 
in  many  applications.  The  most  serious  disadvantage  of  the 
dual  scatter  arrangement  is  the  small  amount  of  collected 
light  and  the  resulting  poor  signal  quality.  Therefore, 
more  powerful  lasers  are  used  to  increase  the  intensity  of 
the  scattered  light.  The  basic  detection  methods  of  the 
Doppler  frequency  for  three  different  modes  are  described  in 
Figures  3,  4 and  5,  respectively.  In  order  to  help  the  under- 
standing of  the  above  Figures,  the  general  relation  between 
scattered  light  waves  and  Doppler  shift  is  shown  in  Figure  2 
[7,  8*].  Among  the  three  basic  modes,  the  dual  beam,  forward 
scattering  system  is  now  most  popular,  since  it  is  generally 
easiest  to  use,  gives  the  best  signal  to  noise  ratio,  and  is 
quite  versatile.  This  is  the  system  that  the  MIT  Marine  Hydro- 
dynamics Laboratory  has  recently  acquired  mainly  for  the 
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study  of  the  flow  around  model  propellers.  Therefore,  the 
working  theory  of  dual  beam,  forward  scattering  system  will 
be  discussed  in  detail. 

(1)  Fundamental  Theory  of  Dual  Beam  Mode 

As  shown  in  Figure  6,  a basic  characteristics  of  the 
dual  beam  mode  is  the  crossing  of  two  laser  beams  of  equal 
intensity  at  a p>int  in  the  fluid  to  be  measured.  To  achive 
this,  the  laser  beam  is  first  separated  into  two  parallel 
beams  of  equal  intensity  that  are  equidistant  from  the 
original  beam.  These  beams  are  then  focused  by  a lens. 

Where  the  two  beams  cross,  they  interfere  each  other  to 
form  alternate  regions  of  high  and  low  intensity  light, 
called  "fringes",  as  shown  in  Figures  6-(c)  and  7.  These 
fringes  are  caused  by  the  light  in  the  two  beams  cancelling 
each  other  in  some  regions  and  complementing  each  other  in 
others.  This  interference  effect  is  described  in  Figure  8. 
As  the  particle  passes  through  the  crossing  point  in  the 
plane  of  two  beams,  these  variations  in  incident  light 
cause  variations  in  the  intensity  of  the  scattered  light 
from  the  particle.  If  this  scattered  light  from  the  parti- 
cle is  picked  up  by  a photodetector,  it  can  be  converted 
to  an  electrical  signal  whose  frequency  is  proportional  to 
the  rate  at  which  the  particle  is  crossing  the  interference 
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fringes.  A typical  photodetector  output  when  observed  on 
an  oscilloscope  is  shown  in  Figure  9.  The  laser  beams 
generally  have  a Gaussion  intensity  distribution.  This 
causes  the  amplitude  variation,  as  shown  in  the  Figure  6-(d), 
with  the  maximum  amplitude  occurring  in  the  center  of  the 
beam  crossing  area.  The  frequency  of  the  signal  in  Figure  9 
is,  of  course,  proportional  to  the  velocity  with  which  the 
particle  is  moving  across  fringes.  As  shown  in  Figure  10, 
the  distance  between  fringes  (fringe  spacing),  df,  is  not 
affected  by  the  refractive  index,  but  depends  only  on  the 
wavelength  of  the  laser  light,  and  the  intersection 
angle  of  the  beams,  2€,  tnrough  the  relation: 


d 


f 


\ 

2 sin€ 


(1) 


Multiplying  this  distance  by  the  output  frequency  (number 
of  fringes  traversed  per  unit  of  time)  will  give  the  parti- 
cle velocity: 


Only  the  intersection  angle  and  laser  wavelength  are  required 
for  interpreting  the  frequency  in  terms  of  velocity.  With 
lasers,  the  wavelength  is  known  accurately  (error  < 0.01%  [9~|) 
and  an  accurate  measurement  of  beam  crossing  angle  is  possible 
for  the  given  optics. 
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Also,  the  only  velocity  measured  is  the  component  perpen- 
dicular to  the  "fringes".  It  should  be  noted  that  a 180° 
ambiguity  exists  in  that  the  system  cannot  distinguish 
"reverse"  flow  from  "forward"  flow.  Technique  for  fre- 
quency shifting  one  beam  to  permit  this  are  discussed  in 
Section  II.3.E.  An  aperture  is  used  on  the  photodetector 
so  only  particles  going  through  the  beam  crossing  area 
will  give  a signal.  Receiving  optics  and  photodetector 
with  the  dual  beam  system  can  be  placed  at  any  angle  and 
still  give  same  frequency.  However,  signal  quality  and 
intensity  will  vary  greatly.  In  practice,  the  receiving 
optics  and  photodetector  is  usually  arranged  in  line  with 
the  focusing  optics,  since  this  is  the  location  of  the 
strongest  signal. 

(2)  Direction  of  Measurement 

The  measurement  direction  is  the  vector  at  right  angles 
to  the  axis  of  the  optical  unit  and  in  the  plane  of  the  two 
beams.  In  general,  any  direction  in  a plane  perpendicular 
to  the  axis  of  the  optical  unit  can  be  selected  simply  by 
rotating  the  beam  splitter.  The  actual  measuring  location 
does  not  move  during  this  rotation.  The  ability  to  rotate 
the  beam  splitter  has  several  applications! 
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1.  It  can  be  rotated  to  determine  the  flow  direction 
(maximum  reading) 

2.  Readings  at  + 45°  to  the  mean  flow  can  be  used  to 
measure  uv  (Reynolds  Stress) 

3.  Readings  at  + 45°  can  be  utilized  to  obtain  an 
accurate  measure  of  the  mean  velocity  vector 

(3)  Volume  of  Measurement 

When  using  LDA,  the  size  of  measuring  volume  affects  the 
accuracy  of  measurements.  It  is  sometimes  very  important  to 
minimize  the  measuring  volume.  Thin  boundary  layers,  small 
test  sections,  accurate  positioning,  and  many  other  experi- 
mental circumstances  make  a "point"  measurement  desirable. 

In  this  regard,  it  is  necessary  to  discuss  the  difference 
between  two  measuring  volumes,  that  is,  geometric  measuring 
volume  and  effective  measuring  volume. 

Geometric  Measuring  Volume 

The  geometric  measuring  volume  is  "shaped"  as  shown 
on  Figure  6-(b).  If  the  diameter  of  the  laser  beam  at  the 
measuring  volume  is  d^  and  the  intersection  angle  is  2€ , 


then  t 
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dl 

<^m  cos  € 


(4) 


For  a diffraction  limited  lens,  one  representation  of  the 
beam  diameter  at  the  focal  point  is  [10,  ll]i 

<d*>G  (5) 

From  Eq.  (5),  a short  focal  length  and  large  original  beam 
diameter  gives  the  minimum  measuring  volume  size. 

Effective  Measuring  Volume 

If  the  initial  laser  beam  is  plane  parallel  and  has  a 
Gaussian  intensity  profile,  the  light  waves  near  the  focal 
waist  are  also  plane  parallel  and  have  a Gaussian  intensity 
profile.  This  wave  and  optical  phenomenon  is  schematically 
shown  in  Figure  11.  Due  to  this  Gaussian  intensity  distri- 
bution, the  shape  of  effective  measuring  volume  is  ellipsoidal 
instead  of  diamond  shape.  This  conclusion  has  been  derived 
based  on  the  signal  at  the  photodetector  rather  than  appearance. 
In  Figure  11,  e-^  quantities  are  used  for  convenience.  They 
are  related  to  the  e~2  quantities  by  [9]: 
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In  equation  (6),  the  factor  of  0.8  is  incorporated  following 
the  recommendation  of  [ll].  This  is  the  expression  of  the 
effective  beam  diameter  at  the  crossing  point  for  a diffrac- 
tion limited  lens.  Therefore,  dimensions  of  the  effective 
measuring  volume  are: 
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The  effective  measuring  volume  is: 


<VOl)E 


” (dt>E 
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(10) 


The  dimensions  (Xm  >£  and  (dm  )g  basically  control  the  measuring 
volume  size.  This,  in  turn,  controls  the  spatial  resolution 
of  the  system. 
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(4)  Number  of  Fringes  in  Measuring  Volume 

Figure  6-(c)  schematically  shows  the  fringes  where  the 
beams  cross.  The  number  of  fringes  is  equal  to  the  number 
of  Doppler  cycles  a single  particle  will  generate.  Since 
dm  is  the  diameter  of  the  measuring  volume,  the  number  of 
fringes  is  given  bys 
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Therefore,  the  geometric  number  of  fringes  is, 
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and  the  effective  number  of  fringes  is. 
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It  is  of  note  that  equations  (12)  and  (13)  indicate  that 
the  number  of  fringes  depends  only  on  the  original  laser 
beam  diameter  and  the  distance  between  splitted  beams. 

Of  course,  equations  (12)  and  (13)  are  valid  only  for  the 
very  center  of  the  beam  crossing  point.  Towards  each  end 
of  the  length,  Zm,  the  number  of  fringes  would  decrease. 
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(5)  Aperture 

An  aperture  is  used  in  front  of  or  as  a part  of  the 
photodetector.  This  aperture  blocks  light  coming  from 
locations  other  than  the  beam  crossing  point  and  thereby 
improves  the  signal  to  noise  ratio  of  the  system.  In 
addition,  it  prevents  the  possibility  of  an  invalid  Doppler 
signal  being  generated  by  mixing  of  light  scattered  by  a 
particle  in  each  of  the  two  beams  but  not  at  the  crossing 
point.  The  aperture  size  is  determined  by  the  magnifica- 
tion which  is  defined  as. 


Therefore,  from  the  previous  equations. 


(14) 


d 


a 


(dmVM 


± F r 
TT  D2  cosf 


(15) 


As  shown  in  equation  (14),  determination  of  the  aperture 
size  depends  on  the  diameter  of  the  beam  exiting  the  laser 
and  the  focal  lengths  of  transmitting  and  receiving  lenses. 
If  the  focal  lengths  of  transmitting  and  receiving  lenses 
are  same,  i .e. , 

Ft  = Fr 

then, 

da  = (dm)E  * (d*)E 

This  is  the  case  of  the  maximum  aperture  efficiency. 
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(6)  Selection  of  Lenses 

The  lens  selection  depends  on: 

1 . Dimension  of  test  section 

2.  Flow  velocity  range 

3.  Frequency  response  of  signal  processor 

The  focal  length  of  both  transmitting  and  collecting  lenses 
must  be  long  enough  to  reach  the  test  point  corrected  for 
windows  and  fluids  other  than  air.  Velocity  range  and 
frequency  response  of  signal  processor  are  the  next  consid- 
eration. Knowing  the  maximum  velocity  and  the  maximum  fre- 
quency response  of  signal  processor,  the  desired  transmit- 
ting lens  can  be  determined  by  the  following  relation: 


(16) 


For  optimum  signal  intensity,  it  is  usually  best  to  use  the 
shortest  focal  length  lens  that  is  consistent  with  the 
measuring  situation  and  signal  processor.  Also,  it  is 
recommended  [9^]  that  transmitting  and  receiving  lenses 
of  the  same  focal  length  be  used  for  most  applications 
in  order  to  obtain  the  maximum  aperture  efficiency,  as 
shown  in  section  II. 3. A. (5). 


34 


B.  Particulates 

The  laser  anemometer  actually  measures  the  velocity 
of  particles  in  the  flow,  not  fluid  velocity,  because  it 
measures  velocity  by  detecting  light  scattered  from  particles 
in  the  fluids  This  fact  is  particularly  important  in  high 
speed,  low  pressure  flows  where  only  small  particles  should 
be  sensed  since  larger  ones  may  not  be  at  the  same  velocity 
as  the  flow.  To  measure  fluid  velocity  with  a laser  anemo- 
meter, therefore,  there  must  be  particles  in  the  flow  in 
sufficient  concentration  and  with  proper  size  to  give  the 
reguired  signal.  The  guality  of  the  data  depends  strongly 
on  the  concentration  and  type  of  particulates  in  the  flow. 

As  well  as  the  number  of  particulates,  uniform  size  particles 
will  also  give  better  results  than  a wide  distribution  in 
size.  This  is  due  to  the  improved  signal  provided  for  the 
signal  processor  and  a reduction  in  the  influence  of  particle 
collosions.  As  a general  rule,  liguids  have  adaquate  con- 
centration and  size  of  particulates  for  essentially 
continuous  measurements.  Also,  the  relatively  high 
density  of  most  liquids  and  the  low  accelerations  minimize 
the  importance  of  particle  size.  For  most  gas  flows,  the 
number  of  particulates  is  insufficient  to  give  data  rates 
high  enough  to  provide  an  essentially  continuous  output. 
Therefore,  it  is  apparent  that  proper  seeding  is  important 
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for  most  measurements  in  gases.  This  is  without  considering 
the  effect  of  particle  size.  It  should  be  noted  that  other 
characteristics  of  the  particles  besides  size  and  concen- 
tration can  influence  the  reading.  These  include  refrac- 
tive index  (important  in  water),  rotation,  collisions, 
density,  and  interaction  with  the  flow.  However,  size 
and  concentration  are  the  dominant  influences  in  terms 
of  the  quality  of  the  data  obtained  with  a laser  anemo- 
meter[l2,  13]. 

C.  Size  and  Type  of  Laser 

Lasers  come  in  a wide  range  of  sizes,  vravelengths  and 
quality.  Tradi tionslly , lower  power  (5  mw  or  15  mw)  He-Ne 
lasers  have  been  used  for  small  laboratory  type  experiment. 
It  has  a long  tube  life,  is  relatively  inexpensive  in  the 
smaller  sizes,  and  is  readely  available  from  a number  of 

sources.  In  the  larger  size  it  is  expensive  and  its  wave- 

« 

length  (6328A)  is  not  optimum  for  most  photomultipliers. 

In  large  tunnels  and  other  difficult  measurement  situations, 
higher  power  (over  100  mw)  Argon-iron  lasers  have  been  used 

O 0 

almost  exclusively.  The  wavelength  (4880A  and  5145A)  are 
in  the  optimum  range  for  most  photomultiplier  tubes.  Other 
lasers  such  as  Krypton,  YAG  and  He-Cd  lasers  might  also  be 
considered  for  some  applications.  For  a specific 
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application,  with  sufficient  data  it  should  bo  possible  to 
compute  laser  power  required  from  the  optical  configuration 
and  characteristics  of  the  fluid,  particulates  and  photo- 
detector. In  practice,  this  is  very  complex  and  can  not 
be  easily  defined  precisely.  However,  some  generalisations 
from  experience  can  be  made  for  initial  selections. 

In  the  dual  beam  forward  scatter  mode,  velocities  below 
200  ft/sec,  adequate  seeding  and  less  than  300  mm  focal 
length  lenses,  very  good  Doppler  signals  are  obtained  with 
good  optics  and  a 5 mw  He-Ne  laser.  Measurements  can  be 
made  in  considerably  more  adverse  conditions,  but  the  signal 
to  noise  ratio  of  the  Doppler  signal  will  be  somewhat  de- 
graded. In  backscatter  for  the  above  conditions,  a larger 
laser  is  required  for  good  Doppler  signals.  A 15  mw  unit 
is  adequate  for  the  shorter  focal  length  lenses.  In  forward 
scatter,  a 15  mw  unit  works  very  well  with  the  600  mm  lenses. 
At  shorter  focal  lengths  in  forward  scatter,  velocities  up 
to  Mach  1 can  be  handled  nicely  with  a 15  mw  laser.  The 
size,  type  and  wavelength  of  commercially  available  lasers 
is  rapidly  increasing  so  the  above  discussions  may  soon  be 
outdated.  To  summarize,  for  most  laboratory  experiments  a 
5 mw  or  15  mw  He-Ne  laser  is  adequate  for  good  measurements. 
This  is  particularly  true  since  forward  scatter  can  generally 
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bo  utilized.  As  backscatter  becomes  necessary,  the  focal 
length  increases  and/or  the  velocity  goes  up,  transition  to 
a larger  laser  such  as  an  Ar-Ion  unit  becomes  necessary. 

All  of  the  above  discussions  are  for  dual  beam  measurements. 
Generally,  reference  beam  measurements  require  somewhat 
more  power  for  equal  signal  to  noise  ratio  while  dual 
scatter  measurements  need  considerably  more  power.  With 
a 5 mw  laser,  dual  scatter  measurements  should  be  attempted 
only  with  short  focal  length,  relatively  low  velocities  and 
forward  scatter.  In  dual  scatter  measurements  in  the 
backscatter  mode,  even  the  15  mw  laser  is  inadequate. 
Therefore,  the  selection  of  laser  depends  not  only  on  the 
size  of  the  laboratory,  but  also  on  the  mode  of  optical 
arrangements . 

D.  Signal  Processors 

The  presently  available  signal  processors  can  be 
devided  into  two  types: 

Special  Techniques 

Spectrum  Analyzer 
Analog  Correlator 
Photon  Count  Correlator 
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Frequency  to  Voltage  Conveters 

Tracker 
Counter 
Filter  Bank 
Direct  to  Computer 

Each  of  the  above  is  the  technique  for  processing  signals 
from  the  photodetector  of  a LDA  system.  The  features  of 
various  techniques  will  be  briefly  discussed. 

( 1 ) Spectrum  Analyzer 

Spectrum  analyzers  were  the  first  instruments  for  ana- 
lyzing LDA  signals.  Spectrum  analyzers  generally  sweep 
through  the  frequency  range  of  interest,  displaying  the 
frequencies  present  as  peak  on  some  type  of  display.  They 
are  readily  available  and  permit  a good  measure  of  the  mean 
velocity  and  an  estimate  of  turbulence  intensity.  Problems 
with  this  method  include* 

1.  data  is  not  "real  time"  (analyzer  must  sweep  through 
bandwidth  of  interest) 

2.  broadening  of  spectrum  due  to  the  finite  filter  bandwidth 

3.  broadening  of  spectrum  due  to  a signal  of  finite  length 


which  includes  high  order  terms 
4.  broadening  of  spectrum  due  to  noise 
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Tho  item  2 above  can  be  corrected.  The  item  3 need  not  be 
a problem  with  most  signal  processors  when  only  one  particle 
at  a time  is  in  the  measuring  volume,  but  will  show  up  on  a 
spectrum  analyzer.  Figure  12  shows  the  basic  idea  of  signal 
processing  method  of  spectrum  analyzer  and  spectrum  broadening. 

(2)  Analog  Correlator 

Correlators  sample  the  signal.  Standard  analog  correlators 
have  been  used  very  little  for  LDA  work,  probably  because  they 
are  too  slow,  particularly  at  high  frequencies.  However, 
within  its  frequency  range,  the  analog  autocorrelator  will 
give  both  mean  velocity  and  a measure  of  turbulence  intensity. 

In  this  regard,  this  information  is  sililar  to  a spectrum 
analyzer . 

(3)  Photon  Count  Correlator 

The  photon  count  correlator  is  basically  the  same  as  the 
analog  correlator  with  a special  input  section  for  pulse 
counting  instead  of  analog  voltage  levels.  It  operates  on 
individual  photon  pulses  and  would  convert  the  pulse  frequency 
to  a digital  word.  The  reason  for  a special  instrument  is  the 
need  for  speed.  The  photon  correlator  will  operate  at  very 
low  light  intensities.  This  is  its  primary  advantage.  It 
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permits  measurements  where  other  techniques  do  not  have  an 
adequate  signal.  If  light  intensity  is  too  high,  individual 
photon  pulses  cannot  be  resolved.  However,  one  can  almost 
always  decrease  light  intensity,  while  the  opposite  is  not 
true.  Mean  flow  and  turbulence  intensity  are  available  from 
this  technique,  but  time  information  on  flow  signal  is  not 
readily  available.  For  frequency  to  voltage  converter  type 
signal  processors,  some  important  characteristics  should  be 
mentioned  prior  to  the  discussions.  They  arei 

1.  number  of  Doppler  cycles  required  for  accurate 
measurement 

2 . capture  bandwidth 

3.  slew  rate 

4.  dynamic  range 

5.  ability  to  extract  signal  from  noise 

6 . measurement  accuracy 

The  items  2 and  4 above  are  particularly  important  from  the 
operational  view  point. 

(4)  Tracker 

Tracker  is  basically  a combination  of  tracking  filter 


and  frequency  to  voltage  converter.  Typically,  this  type  of 
signal  processors  has  been  used  where  the  particle  concen- 
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, tration  is  hiqh,  since  it  is  easier  to  use  with  multiple 

t particles  in  the  measuring  volume  and  with  relatively  noisy 

signals.  The  fundamental  limitation  of  any  tracker  type 
signal  processors  is  the  "Capture"  range.  If  the  change  of 
velocity  between  one  particle  and  the  next  particle  is 
greater  than  the  "Capture"  range  of  the  unit,  it  will  not 
track  the  new  value  and  any  succeeding  values  that  remain 
outside  the  capture  range.  This  happens  particularly  when 
a combination  of  low  particle  concentrations  and  large 
changes  in  velocity  occur.  Figure  13  illustrates  the  dynam- 
ic range  and  capture  bandwidth.  From  Figure  13  the  follow- 
ing operational  points  can  be  derived: 

1.  operation  should  be  towards  full  scale  for  optimum 
extraction  of  the  signal  from  noise 

2 . operation  should  be  towards  the  bottom  of  the  range 
for  maximum  capture  bandwidth 

With  high  particle  concentration  in  the  fluid,  noise  due  to 
varying  particle  phase  will  be  present.  A low  pass  filter 
should  always  be  used  to  minimize  this  noise  by  filtering  out 
frequencies  above  those  of  interest.  Unfortunately,  all  of 
this  phase  noise  can  not  be  filtered  out,  since  it  is  essen- 
tially "white"  noise.  Practically,  there  are  limits  to  both 
accuracy  and  minimum  output  noise  levels  on  trackers  even  with 
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ideal  input  signals,  and  this  greatly  restricts  the  turbu- 
lence measurements  when  the  intensity  is  high. 

( 5)  Counter 

A counter  measures  the  time  interval  for  a fixed  number  of 
cycles,  N,  of  the  Doppler  signal.  The  inverse  of  this  time 
interval  is  then  velocity.  This  type  of  signal  processors 
has  been  applied  in  the  flow  situations  of  high  velocity, 
low  particle  concentration.  It  is  often  used  in  difficult 
measuring  situations  with  the  large  laser.  A counter  does 
not  have  a "Capture"range , since  any  signal  within  its  entire 
freguency  range  can  be  "Captured".  The  accuracy  of  the  system 
mostly  depends  on  the  accurate  zero  crossing  measurements. 

For  a single  particle  and  low  background  noise,  N cycle  of 
the  Doppler  signal  represent  the  correct  frequency  for  the 
velocity  present.  With  multiple  particles,  phase  changes 
can  give  a significant  error  over  N cycle,  particularly  if 
N is  small.  Also,  the  "pedestal"  or  DC  component  of  the 
Doppler  burst  must  be  removed  for  accurate  zero  cross  count- 
ing. Ideally,  it  can  be  done  without  limiting  the  dynamic 
range  of  the  instrument.  Generally,  the  following  methods 
are  utilized  to  remove  the  pedestal j 
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1.  high  pass  filter 

2.  optical  pedestal  removal 

3.  frequency  shifting 

4.  filter  bank 

If  the  method  of  filtering  is  used  to  removed  the  pedestal, 
the  frequency  range  and  hence  the  "Capture”  range  will  be 
limited  by  the  filter.  Furthermore,  improper  filtering  may 
be  the  another  source  of  error.  Figure  14  schematically 
describes  the  effect  of  noise  on  zero  crossing  measurements. 

A comparison  of  some  characteristics  between  presently  avail- 
able commercial  tracker  type  and  counter  type  signal  proces- 
sors is  shown  in  Table  1.  This  comparison  has  been  made 
because  they  are  the  most  common  techniques. 

(6)  Filter  Bank 

The  filter  bank  consists  of  many  individual  filters. 
Output  essentially  looks  for  filter  that  is  resonating  and 
identifies  it  as  frequency.  The  problems  in  this  type  of 
signal  processors  arei 

1.  several  cycles  are  required  for  resonance 

2.  accuracy  is  limited  by  the  number  of  filters 

The  filter  bank  signal  processors  are  effectively  real  time 
spectrum  analyzers. 

(7)  Direct  to  Computer 


A fast  analog  to  digital  converter  with  subsequent  computer 


data  processing  can  also  bo  used.  The  converter  essentially 
puts  the  "raw"  Doppler  signal  into  digital  form  for  use  in  a 
computer.  The  data  is  then  sorted  out  by  proper  programming. 
Disadvantages  of  this  kind  of  signal  processing  method  are: 

1.  limited  frequency  of  analog  to  digital  converter 

2.  large  amounts  of  computer  time  and  data  storage  required 

Summary 

If  signal  strength  is  adequate,  most  well  designed  signal 
processors  will  work  for  most  flows.  It  is  the  author's 
opinion  that  the  special  techniques  work  better  in  turbulence 
measurements,  although  they  do  not  give  an  instantaneous 
velocity  value.  Furthermore,  photon  correlators  will  operate 
at  very  low  light  levels.  The  tracker,  counter  and  filter 
bank  are  all  methods  to  obtain  a nearly  "instantaneous" 
conversion  of  the  frequency  to  voltage.  Tracker  can  give 
better  results  for  lower  signal/noise  ratio,  are  easier  to 
use  and  cost  less  than  counter  type  processors.  In  highly 
turbulent  flows  with  low  particle  concentration  (mostly 
gases  with  natural  seeding),  a tracker  may  not  work  at  all 
and  a counter  must  be  used.  Counters  will  handle  a wider 
variety  of  flows  than  trackers.  The  above  comments  are 
related  to  the  signal  processors  now  available.  Certainly, 
the  practical  implementation  of  both  techniques  will  be 
greatly  improved  in  the  future. 
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E.  Data  Analysis 

The  fact  that  laser  anemometer  actually  measures  the 
particle  velocity,  not  fluid  velocity,  gives  special  problems 
in  data  analysis.  Particles  are  distributed  randomly  in  the 
fluid  and  the  naturally  occurring  particles  vary  in  size  and 
concentration.  Intentional  seeding  permits  control  of  size 
and  concentration  but  still  leaves  a random  spatial  distribu- 
tion. Perhaps  this  is  the  single  most  important  complication 
of  the  LDA  technique.  In  addition,  the  signal,  as  far  as  any 
realizable  signal  processor  is  concerned,  is  essentially 
discontinuous  even  with  a high  particle  concentration  due  to 
widely  varying  amplitude.  In  many  applications,  the  signal 
is  actually  intermittent  due  to  the  absence  of  particles  in 
the  measuring  volume.  The  choice  of  the  best  method  of  data 
analysis  depends  on  many,  often  interrelated  factors  such  asj 

Particle  concentration  (Burst  Density) 

Particle  size 

Spatial  / Temporal  Resolution 

Signal  to  noise  ratio 

Signal  processing  method 

/ 

Signal  validation  method 

Type  of  information  required 


Among  the  factors,  the  single  most  important  one  is  the 
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particle  concentration.  The  statistics  of  high  density 
signals  are  fundamentally  different  from  low  density  sig- 
nals. Most  data  analysis  systems  such  as  rms  meters, 
spectrum  analyzers,  correlators,  and  even  voltmeters  for 
mean  flow  are  designed  to  operate  on  continuous  data. 

Therefore,  it  is  ideal  when  the  output  of  the  singal  pro- 
cessor on  a laser  anemometer  can  be  considered  continuous. 

This  depends  primarily  on  the  data  rate.  For  mean  flow, 
the  fact  that  the  processor  holds  the  last  reading  helps 
considerably.  When  this  happens,  there  are  at  least  two 
ways  of  time  averaging  the  data.  One  way  is  to  average  the 
signal  over  all  time  regardless  of  the  signal  validity.  This 
corresponds  to  low  pass  filtering  the  data.  The  other  way  is 
to  average  the  signal  over  only  those  intervals  in  which  the 
signal  was  valid.  While  this  approach  is  intuitively  more 
correct,  it  is  also  more  difficult  to  implement.  In  the  case 
of  high  burst  density  signals,  there  is  little  difference  be- 
tween the  two  types  of  averages.  For  most  purposes,  therefore, 
it  is  sufficient  to  approximate  the  latter  by  the  former 
for  high  burst  density  signals.  However,  this  may  not  be 
true  of  low  burst  density  signals.  For  turbulent  flow, 
true  rms  is  normally  the  first  parameter  measured  after  the 
mean  flow.  Here,  the  validity  of  assuming  a continuous  meas- 
urement depends  on  the  data  rate  and  the  frequencies  that 
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contain  the  primary  energy.  Certainly,  on  the  average, 
several  data  points  per  cycle  are  required  for  frequencies 
that  represent  a significant  part  of  the  total  turbulent 
energy.  For  most  turbulence,  most  of  the  energy  is  in  the 
lower  frequencies  so  that  mean  velocity  and  turbulence 
intensity  can  frequently  be  measured  assuming  high  data 
rates.  The  use  of  continuous  data  procedures  for  spectral 
analysis  depends  directly  on  the  number  of  data  points. 
Again,  at  the  maximum  frequency  of  interest,  it  is  necessary 
to  have,  on  the  average,  several  data  points  per  cycle  for 
good  data.  For  correlation  with  time  dalay,  the  correlation 
curve  will  be  biased  if  a new  data  point  has  not  occurred  at 
the  minimum  time  increment  utilized  in  the  correlation. 

A detailed  analysis  including  both  the  samples  per  second 
and  the  distribution  of  these  samples  is  necessary  to  pre- 
cisely define  the  validity  of  the  data  when  assuming  a con- 
tinuous output.  The  fact  the  sample  points  are  not  equally 
spaced  can  actually  be  utilized  to  improve  the  effective 
response  without  seperately  keeping  track  of  time  intervals 
if  an  appropriate  estimator  is  used[l4].  For  low  burst 
density  signals,  the  best  method  of  taking  data  is  to  take 
individual  data  points  as  they  occur.  It  is  important  to 
notice  that  low  burst  density  means  a low  probability  of 
more  than  one  particle  in  the  measuring  volume  at  a time. 


48 


i 

i 


There  still  can  be  a sufficiently  high  data  rate  to  follow 
the  flow  characteristics  of  interest.  At  low  data  rates, 
maximum  data  utilization  is  obtained  when  both  time  data 
and  velocity  data  are  utilized.  Here,  the  term,  "low  data 
rate”,  generally  means  a rate  so  low  that  even  the  large 
scale  fluctuations  in  flow  are  not  recognizable  on,  say, 
an  oscilloscope  trace  of  the  signal  processor  output. 

With  an  on-line  computer,  the  computer  clock  can  be  used 
to  measure  the  time  interval  between  data  points.  Combining 
this  with  the  velocity  measurement  represented  by  each  data 
point  permits  complete  statistical  analysis  of  the  results. 

F.  Other  Systems 

There  are  numerous  modifications  of  the  three  basic 
systems  that  can  be  utilized  for  special  measurement 
requirments.  These  may  indicate  the  flexibility  of  the 
laser  anemometer.  However,  only  the  frequency  shifter 
will  be  discussed  here. 

Frequency  Shifting 


The  standard  laser  anemometer  system  has  a 180°  direction 
ambiguity.  In  other  words,  it  cannot  tell  whether  the  flow  is 
forward  or  backward.  This  becomes  a limitation  in  a number  of 
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measurements.  In  order  to  expand  the  applications  of  laser 
anemometry,  a frequency  shifter  is  used.  Proper  application  of  the 
frequency  shifter  permits  measurements  in  reverse  flow  as  well 
as  assistinq  in  other  difficult  applications.  In  principle, 
it  can  be  utilized  with  any  mode.  In  practice,  it  has  been 
used  primarily  with  the  dual  beam  and  reference  beam  modes. 

In  the  reference  beam  mode,  the  frequency  of  the  scattered 
light  is  compared  with  that  of  a reference  beam.  At  zero  flow, 
the  two  light  frequencies  are  equal  and  the  frequency  output 
of  the  photodetector  (difference  between  reference  beam  fre- 
quency and  scattered  light  frequency)  goes  to  zero.  If  the 
reference  beam  is  frequency  shifted,  for  example,  40  MHZ, 
then  the  photodetector  output  at  zero  flow  is  40  MHZ.  The 
frequency  output  will  be  higher  for  particle  velocities  in 
one  direction  and  lower  for  the  other  direction.  Essentially, 
the  resulting  system  gives  a zero  offset  that  permits  bi- 
directional measurements  with  a laser  anemometer.  In  the 
dual  beam  mode,  two  beams  of  the  same  frequency  interfere  with 
each  other  where  they  cross  and  this  interference  causes  a 
stationary  fringe  pattern.  If  one  beam  is  shifted  in  fre- 
quency relative  to  the  other,  the  fringes  can  be  thought  of 
as  "moving".  If  the  difference  in  frequency  of  the  two  beams 
is  40  MHZ , a stationary  particle  in  the  measuring  volume  will 
then  give  a 40  MHZ  signal  from  the  photodetector  picking  up 
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the  scattering  light.  Movement  in  the  same  direction  as  the 
moving  fringes  will  lower  the  frequency  while  movement  in  the 
opposite  direction  will  raise  it.  Again,  the  result  is  equiva- 
lent to  a zero  offset.  In  short,  this  technique  can  be  sum- 
marized as  "shifting  the  frequency  of  the  light  of  one  of  the 
heterodyne  signals  by  fs  so  that  zero  velocity  will  give  fre- 
quency fs  at  the  photodetector.  Techniques  for  frequency 
shifting  include  acousto-optic  (Bragg)  cells  and  rotating 
diffraction  gratings.  Both  techniques  work  well,  but  acousto- 
optic (Bragg)  cell  technique  is  more  widely  used  due  to  the 
small  ranqe  of  maximum  shift  frequency  of  present  rotatinq 
diffraction  grating  technique.  The  principles  of  the  Debye- 
Sears  effect  form  the  basis  of  acousto-optics  cells,  which  may  be 
used  to  divide  a light  beam  into  several  beams  with  different 
frequencies  and  directions.  Bragg  cell  is  a particular  case 
of  these  cells  operating  to  yield  one  beam  with  a frequency 
and  direction  different  from  the  incident  beam.  As  shown  in 
Figure  15 , the  incident  light  beam  passes  through  a medium 
in  which  acoustic  waves  are  travelling.  Depending  on  the 
angle,  the  acoustic  frequency  (normally  40  MHZ)  is  either 
added  or  subtracted  from  the  light  frequency.  The  intensity 
of  the  light  in  each  refracted  beam  depends  on  several  para- 
meters including  the  cell  dimensions,  the  acoustic  and  light 
frequencies,  the  absorbing  medium,  the  electrical  power  used 
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to  drive  the  crystal  and  the  angle  between  the  acoustic  wave 
fronts  and  the  incident  beam.  In  practice,  the  intensity  of 
the  two  beams  coming  out  (zeroth  order  and  first  order)  can 
be  equalized.  Frequency  shifting  is  a very  useful  tool  in 
LDA.  The  specific  applications  of  frequency  shifting  are: 

1.  Flow  reversals  in  the  mean  flow  direction 

2.  Turbulence  components  perpendicular  to  the  mean  flow 

3.  Pedestal  removal  by  frequency  shifting 

4.  Optimizing  the  frequency  range  on  signal  processor 

II. 4.  Errors  in  Laser  Doppler  Anemometry 

t 

The  errors  that  are  associated  with  measurements  using 
laser  anemometry  come  from  many  theoretical  and  practical 
limitations.  They  can  be  listed  as: 

A.  velocity  fluctuations  within  the  volume 

' > 

B.  mean  velocity  gradients  across  the  volyme 

C.  finite  transit  time  of  the  particles  through  the  volume 

D.  low  particle  concentration 

E.  signal/noise  ratio 

F.  precision  of  the  signal  processors 

The  items  A,  B and  C above  are  theoretical  limitations  associ- 
ated  with  finite  measuring  volume  [15,  16].  The  items  A and  B 
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effects  similar  to  those  found  in  using  a hot-wire  anemometer. 

A.  Fluctuations  in  particle  movement  within  the  measuring 
volume  can  influence  the  resulting  mean  and  turbulence 
data . 

B.  Likewise,  a finite  measuring  volume  can  affect  mean 
flow  measurements  in  steep  velocity  gradients. 

Furthermore,  for  uniform  particle  concentrations,  the 
number  of  data  points  will  be  higher  at  high  velocities. 
Summing  individual  data  points  will  lead  the  result  to 
the  high  velocity  side. 

C.  The  effect  of  finite  transit  time  of  the  particles 
through  the  volume  basically  derives  from  the  fact 
that  when  a spectrum  analysis  is  done  on  the  output 
of  the  photodetector  of  an  LDA,  there  is  a spectral 
broadening.  If  one  considers  that  the  Doppler  signal 
for  a single  particle  has  a beginning  and  an  ending, 
this  in  itself  must  generate  higher  frequency  variations 
that  show  up  as  spectral  broadening  even  in  laminar  flow. 
This,  then  gives  a "minimum''  turbulence  intensity  that 
can  be  measured.  How  much  spectral  broadening  affects 
the  various  types  of  signal  processors  is  somewhat 
controversial.  With  a very  noise  free  Doppler  signal. 
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a system  measuring  time  between  zero  crossing  (counter 
type)  would  seem  to  be  independent  of  the  transit  time 
ambiguity  when  only  one  particle  is  in  the  control 
volume  at  any  one  time.  Similarly,  with  a similar 
signal,  a system  reading  the  signal  while  it  still  has 
good  amplitude  and  holding  that  reading  (tracker  type) 
would  also  seem  to  be  free  of  the  transit  time  ambiguity. 
However,  there  is  probably  never  only  one  particle  in 
the  control  volume  even  though  one  dominates  due  to  size. 
Also,  one  never  has  a truly  noise  free  signal. 

The  items  D,  E and  F are  practical  limitations. 

D.  Low  particle  concentration  will  affect  the  time  averaged 
mean  velocity.  A simple  average  of  output  data  will  not 
give  the  correct  average  value. 

E.  The  effect  of  poor  signal/noise  ratio  is  the  most  impor- 
tant practical  limitation  on  measurement  accuracy.  The 
signal  to  noise  ratio  of  the  Doppler  signal  affects  the 
accuracy  of  any  signal  processor,  particularly  for  measure- 
ments of  turbulent  flow.  In  the  counter,  the  effect  of 
noise  is  clearly  seen  when  either  the  data  rate  goes 

down  or  the  number  of  count  is  decreased.  In  the  tracker, 
a poor  signal  noise  ratio  also  results  in  a decreased 
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sample  rate  as  well  as  some  degradation  in  the  track- 
ing accuracy  and  stability.  There  are  numerous  sources 
of  noise.  They  can  be  divided  into  two  categories: 

Inherent  Sources  of  Noise 

(1)  Variations  of  fluid  velocity  across  the 
measuring  volume 

(2)  Multiple  particle  signals 

(3)  Shot  noise 

(4)  Brownian  motion 

Extraneous  Sources  of  Noise 

(5)  Reflected  light  waves  from  the  laser 

(6)  Light  waves  from  the  environment 

(7)  Noise  from  frequency  shifting 

(8)  Optical  noise 

(9)  Signal  processor  noise 

(1)  Noise  from  the  variation  of  fluid  velocity 

This  is  a source  of  noise  associated  with  the  measuring 
volume  size.  This  is  one  of  some  unique  sources  of  the  laser 


anemometer,  when  compared  with  the  hot  wire  anemometer.  This 
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noise  source  is  not  present  with  a hot  wire  which  averages 
variations  of  heat  transfer  over  its  length.  Whereas  the 
laser,  measuring  on  individual  particles,  reacts  much 
differently . 

(2)  Noise  from  multiple  particle  signals 

This  is  a second  inherent  source  of  noise  - commonly 
called  "Doppler  ambiguity"[l6] . The  combination  of  varying 
Doppler  signal  amplitudes  and  random  phases  between  signals 
from  different  particles  results  in  artificial  variations 
in  the  freguency  seen  by  the  signal  processor.  This  noise 
is  not  a problem  at  low  particle  concentrations,  but  becomes 
significant  when  the  probability  of  multiple  particles  in 
the  measuring  volume  is  high.  This  phenomenon  particularly 
limits  the  minimum  turbulence  that  can  be  measured.  Several 
methods  are  available  to  reduce  the  effect  of  Doppler  ambi- 
guity. A low  pass  filter  should  be  used  so  that  only  the 
freguencies  of  interest  are  observed.  As  discussed  previously, 
» however,  all  of  this  noise  cannot  be  filtered  out  since  it  is 

essentially  "white"  noise.  Cross  correlating  two  independent 
measurements  at  the  same  point  in  the  fluid  will  give  a cor- 
relation independent  of  the  ambiguity  noise  [17].  From  this, 
a transform  can  be  used  to  obtain  spectrum.  In  addition  to 
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filtering,  an  alternate  technique  has  been  suggested[ 18] . 

Since  the  ambiguity  spectrum  is  equivalent  to  white  noise, 
its  influence  can  be  determined  by  measuring  the  rms  of  the 
signal  using  different  low  pass  positions.  Once  established, 
it  can  be  subtracted.  Among  the  three  methods  discussed, 
the  most  common  and  popular  method  is  the  use  of  a filter. 

(3)  Shot  Noise 

This  is  the  noise  associated  with  the  discrete  emission 
of  photo-electrons.  If  the  only  light  hitting  the  photo- 
detector is  scattered  by  a single  signal  producing  particle, 
the  dark  current  in  the  photodetector  is  zero  and  there  is 
no  noise  in  the  electronics  following  the  photodetector,  then, 
the  only  remaining  source  of  noise  is  the  shot  noise  associ- 
ated with  the  electron  current  flow  produced  by  the  signal 
itself.  The  mean  square  of  the  shot  noise  fluctuations  is 
proportional  to  the  mean  photo  current.  This  signal-to-shot 
noise  ratio  serves  as  a measure  of  the  performance  of  a LDA, 
and  it  is  the  first  priority  to  insure  adequately  large  values 
of  this  in  designing  or  specifying  a LDA. 

(4)  Noise  due  to  Brownian  motion 

This  can  be  another  possible  source  of  noise.  As  briefly 
discussed  in  reference  6,  this  effect  can  be  a problem  when 
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flow  velocity  is  very  low.  The  molecular  speed  due  to 
Brownian  motion  is  proportional  to  the  following  relation, 

V ~ 7 T/M  (17) 

For  gases,  molecular  weight  is  so  small  that  the  molecular 
speed  is  enormous  and,  hence,  the  frequency  at  the  measuring 
volume  due  to  this  motion  often  reaches  the  order  of  ICh*  MHz. 
This  is  the  reason  why  LDA  measurement  is  not  possible  with- 
out particles.  For  liquids,  this  effect  is  much  less  than 
for  gases,  but  still  gives  high  order  of  frequency. 

(5)  and  (6)  Noise  due  to  reflected  light  waves  from  the  laser 
and  light  waves  from  the  environment 

These  items  are  related  to  the  components  of  light  waves 
at  the  photo  detector.  The  total  electric  vector  at  the  photo 
detector  may  consider  oft 

a)  light  waves  scattered  from  one  or  more  illuminating 
beams  by  one  or  more  particles 

b)  "reference"  light  waves  from  the  laser 

c)  extraneous  light  waves  from  the  laser  reflections 

d)  extraneous  light  waves  from  the  environment  (i.e.  room 
light,  radiation  from  flames,  etc.). 

Each  type  of  light  wave  contributes  to  the  photo  current,  but 
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only  a)  and  b)  produce  signals  containing  useful  velocity 
information.  Laser  light  reflected  or  scattered  by  walls  may 
interfere  with  light  scattered  by  particles  to  produce  noise. 
This  effect  is  particularly  severe  when  measurements  are  made 
very  close  to  walls.  Also,  light  waves  from  the  environment 
can  interfere  with  the  signal  producing  light  waves.  It 
occurs  often  in  practice  particularly  when  the  intensity  of 
scattered  light  waves  by  particles  is  weak. 

(7)  Noise  from  freguency  shifting 

As  discussed  previously,  acousto-optic  (Bragg)  cell 
technique  is  more  widely  used  among  two  common  techniques 
for  frequency  shifting.  However,  Bragg  cell  generates  some 
noise  due  to  high  frequency  vibration  (40  MHZ)  of  local 
oscillator  to  produce  travelling  acoustic  waves.  Usually, 
this  contribution  to  overall  noise  effect  is  regarded  as 
small . 

(8)  Optical  noise 

Poor  optics  design  can  be  another  source  of  noise.  The 
poor  quality  optics,  for  example,  can  generate  the  unparallel 
fringes  or  cause  the  beams  to  cross  outside  the  working  range 
in  the  dual  beam  system.  If  the  fringes  are  not  parallel,  then 
a particle  going  through  different  portions  of  the  measuring 
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volume  will  give  different  frequencies  even  when  traveling 
the  same  velocity.  The  result  is  the  same  as  a mean  velocity 
gradient.  In  addition  to  the  above  effects,  the  optics  whose 
surface  is  not  sufficiently  smooth  would  cause  the  same  noises 
at  those  discussed  in  items  (5)  and  (6). 

t 

(9)  Signal  processor  noise 

The  last  noise  source  to  be  discussed  here  is  the  signal 
processor  itself.  The  minimum  possible  noise  is  that  with  a 
continuous  signal  from  a generator  as  the  source.  This  gives 
a measure  of  the  minimum  attainable  output  noise  with  a given 
signal  processor.  Often,  a discontinuous  signal  that  is  ampli- 
tude modulated  will  give  more  output  noise  than  a continuous 
signal.  This  is  a good  representation  of  the  Doppler  signal 
when  the  particle  concentration  is  light.  A heavy  particle 
concentration,  including  Doppler  ambiguity,  is  difficult  to 
simulate  with  signal  generators.  However,  in  this  case  the 
signal  processor  is  not  usually  the  major  source  of  the  noise. 

The  processor  can  reduce  the  noise  through  either  not  following 
the  rapid  frequency  changes  or  through  output  low  pass  filtering. 

F.  Sometimes,  the  precision  of  the  signal  processors  can  be 
an  important  practical  limitations.  For  example,  the  re- 
lation between  input  frequency  and  output  voltage  should 
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be  linear  for  tracker  or  counter  type  siqnal  processors. 

Since  this  effect  is  so  obvious  that  no  further  discussion 

will  be  necessary. 

S umma  ry 

For  turbulence  or  RMS  measurement,  all  of  the  above  limi- 
tations are  all  interrelated  and  precise  quantitative  error 
analysis  is  practically  impossible.  For  mean  velocity  meas- 
urement, the  error  due  to  phase  noise  may  be  systematically 
reduced  for  periodic  signals  by  averaging  signals  of  successive 
periods  (waveform  analysis).  The  true  signals  add  while  the 
noise  reduces  to  zero  on  average.  In  most  cases  of  mean  velo- 
city measurement  with  reasonable  signal/noise  ratio,  therefore, 
the  following  items  are  considered  to  be  major  factors  which 
affect  the  accuracy: 

1 . accuracy  of  laser  wavelength 

2.  accuracy  of  beam  crossing  angle 

3.  accuracy  of  frequency  shifting  system 

4.  precision  of  signal  processors 

For  most  presently  available  commercial  systems,  the  usual 
error  ranges  are  0.01%,  0.1%.  0.01%  and  0.4%  for  the  above 
items  1,  2,  3 and  4,  respecti vely[ 19] . from  the  above  in- 
formation, therefore,  the  accuracy  of  most  mean  velocity 
measurement  is  expected  to  be  within  99%. 


61 


II. 5.  MIT  Variable  Pressure  Water  Tunnel  and  LDA  Systems 

In  spite  of  the  remarkable  potentials  and  advantages, 
however,  the  actual  application  of  the  Laser-Doppler  Anemo- 
motry  technique  to  the  propeller  tunnel  is  greatly  restricted 
by  the  characteristics  of  tunnel  itself.  Fortunatbly,  the  MIT 
water  tunnel  is  extremely  ideal  for  laser  measurements  for 
the  following  reasons i 

1.  test  section  is  square  with  complete  accessibility 
on  all  four  sides. 

2.  the  size  of  the  test  section  is  large  enough  to 
accomodate  hydrodynamic  models,  and  small  enough 
to  allow  the  use  of  low-powered  lasers. 

3.  the  optical  systems  can  be  rigidly  mounted  on  a 
single  movable  base  which  passes  underneath  the 
test  section. 

MIT  water  tunnel  has  the  following  characteristics j 

Section  shape  j Square  Section 

Inside  dimension  t 20" 

Outside  dimension  « 24" 

Tunnel  window  , i 44"  x 16"  x 2"  plexiglass 

Free  stream  velocity  range  i 0 - 30  FPS 
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The  MIT  water  tunnel  LDA  system  including  signal  process 
and  data  analysis  equipment  has  the  capability  of  meas- 
uring and  storing  the  velocity  at  a fixed  point  in  space 
at  a large  number  of  time  intervals  between  two  blade 
passages.  Combining  these  results  for  a sequence  of 
fixed  measurement  points,  it  is  possible  to  construct 
the  velocity  field  in  propeller  fixed  coordinates  as 
shown  in  Figure  16.  The  MIT  water  tunnel  LDA  system  and 
some  of  its  characteristics  are  shown  in  Tables  2 and  3, 
and  Figures  17  and  18.  In  addition  to  the  tunnel 
characteristics,  there  is  another  important  point  in 
the  practical  application  of  LDA  system.  It  is  the 
quality  of  tunnel  window  (particularly  transmitting  side 
window)  in  the  sense  of  thickness  distribution  and  surface 
smoothness.  As  far  as  thickness  problem  is  concerned, 
what  is  important  is  not  the  absolute  thickness,  but  the 
uniformity  of  thickness.  The  thickness  of  window  should 
be  sufficiently  uniform  for  the  accurate  positioning  of 
beam  crossing  point  and  for  the  perfect  crossing  of  the 
beams.  Furthermore,  its  surface  should  be  sufficiently 
smooth  to  prevent  excessive  surface  scattering  of  the 
laser  beams.  Sample  measurements  of  tunnel  window  thick- 
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ness  distributions  are  shown  in  Figure  19.  Figure  20 
shows  the  relations  between  distances  for  two  different 
transmitting  lenses  and  comparisons  between  some  measured 
and  calculated  distances. 

II. 6.  Calibration 

In  section  II. 4.,  it  has  been  concluded  that  the 
overall  error  of  the  system  would  not  be  greater  than 
1%.  To  check  the  validity  of  this  conclusion,  two  com- 
parisons were  made.  The  first  one  was  the  comparison 
of  tunnel  free  stream  measurement  between  LDA  system 
and  existing  device  (manometer) , and  some  typical  results 
are  shown  in  Table  4.  Since  this  is  the  simplest  case  of 
flow  measurement  situations,  both  should  work  well.  As 
shown  in  Table  4,  the  agreement  between  two  different 
measurements  is  excellent.  In  order  to  make  a comparison 
for  more  complicated  measurements,  a special  device, 
so-called  "calibration  rotor"  was  prepared.  This  rotor 
was  installed  on  the  propeller  shaft  and  rotated  at 
several  different  RPM  to  measure  the  rotating  speed. 
Dimension  of  the  rotor  arrangement  and  results  are 
shown  in  Table  5.  As  shown  in  Table  5,  the  agreement  is 
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very  good  at  lower  RPM  (lower  than  1200).  At  higher 
RPM  (higher  than  1200),  however,  agreement  is  poor  and 
the  measurements  by  two  different  methods  are  contin- 
uously diverging  with  increasing  RPM.  It  was  really 
embarrassing  and  hard  to  explain  at  first.  However, 

it  was  found  that  the  RPM  meter  (counter)  did  not  work  \ 

properly  above  1200  RPM.  This  explained  everything  and 
the  accuracy  of  LDA  measurement  was  confirmed. 

III.  Propeller  Model 

A set  of  four  NSRDC  research  propellers  were  chosen 
for  this  study.  They  are  Propellers  4381,  4382,  4383  and 
4498,  and  were  designed  by  lifting-surface  methods  with 
varying  skew  and  warp[^20,  21  ].  These  propellers  are  con- 
sidered to  be  the  first  model  marine  propellers  designed 
to  study  the  effects  of  skew  systematically.  The  first 
three  propellers  belong  to  one  series,  which  were  designed 
for  uniform  flow  with  all  parameters  held  constant  except 
skew(and  pitch  and  camber  corrections  due  to  skew)[20]]. 
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The  propeller  4498  belongs  to  the  other  series[2l],  which 
differs  from  the  previous  series  in  the  view  point  of  rake. 
However,  all  of  these  propellers  have  the  same  radial  dis- 
tribution of  chord  length  and  thickness,  and  were  designed 
to  have  the  same  loading  distribution  for  uniform  flows 
at  the  design  advance  condition  of  J = 0.889.  Therefore, 
this  set  of  propellers  is  ideal  in  the  study  of  the  effect 
of  large  changes  in  skew  and  rake  without  extraneous  in- 
fluences. The  detail  geometric  and  performance  character- 
istics of  these  propellers  can  be  found  in  references  20 
and  21.  For  convenience,  however,  the  principal  design 
and  some  performance  characteristics  are  reproduced  in 
Tables  4 and  5. 

IV.  Results  of  Experiment 
IV. 1.  General 

The  earliest  and  simplest  measurements  were  made  early 
in  December,  1976,  when  all  of  the  necessary  components  had 
arrived  and  were  assembled  on  a temporary  wooden  base.  These  were 
measurements  of  tunnel  free  stream  velocities.  Since  then, 
a successful  procedure  of  measuring  technique  has  been 
developed  and  established  through  a series  of  trial  meas- 
urements such  as  rotating  speed  of  a specially  prepared 
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rotor,  tunnel  wall  boundary  layer  and  some  field  point 
velocities  for  several  different  propellers.  The  pre- 
liminary test  results  had  not  only  showed  the  accuracy 
and  excellent  repeatability,  which  may  be  the  two  most 
important  factors  in  an  experiment,  but  also  confirmed  the 
possibility  of  successful  measurements  of  the  precise 
quantitative  flow  informations  around  propeller  blades, 
which  would  be  extremely  difficult  or  unsuccessful  with 
conventional  measuring  devices.  The  preliminary  results 
are  partly  reported  in  reference  19  and  will  not  be  dis- 
cussed here.  Only  the  information  directly  obtained  or 
derived  from  a series  of  laser  measurements  for  the  pro- 
pellers introduced  in  Chapter  III  will  be  discussed  in 
the  subsequent  sections.  In  this  study,  results  are 
presented  following  the  coordinate  system  adopted  at 
1973  ITTC  convention,  which  is  shown  in  Figure  21. 

IV. 2.  Slipstream  Radius 

The  conventional  way  of  measuring  slipstream  radius 
basically  relies  on  flow  visualization.  In  this  method, 

propeller  tip  vortices  are  first  made  visible  by  a suitable 

reduction  in  tunnel  pressure  and  radius(or  diameter)  itself 

is  next  measured  either  by  a sighting  telescope  focused  on 
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the  tip  vortices  or  by  photographing  them.  This  visual 
method  is  regarded  as  a reasonable  one,  but  not  the  best 
method,  since  it  has  been  noticed  that  the  shape  of  tip 
vortices  are  somewhat  affected  by  the  tunnel  pressure. 

With  the  LDA  system,  however,  this  can  bo  directly  measured 
without  reducing  the  tunnel  pressure.  In  this  measurement, 
it  has  been  found  that  the  f luctuational  transition  from 
free  stream  to  slipstream  is  extremely  sharp.  This  transi- 
tion occurs  within  the  width  of  only  about  two  percent  of 
propeller  radius  depending  on  the  loading  condition.  Since 
this  information  forms  a basis  of  several  other  derived 
quantities,  the  measurements  were  carefully  performed  at 
two  or  three  different  positions,  for  instance,  approaching 
from  the  top  (along  the  negative  z-axis  with  y = 0)  and  from 
the  side  (along  the  + y-axis  with  z = 0) , and  averaged. 

For  most  cases,  however,  they  are  surprisingly  close  each 
other  and  averaging  was  not  really  necessary.  The  slipstream 
radius  measurements  are  shown  in  Table  8 and  Figures  22 
through  25.  It  may  be  noted  from  the  above  Figures  that 
most  of  contraction  occurs  within  the  axial  distance  of  one 
radius  from  the  propeller  tip.  More  specifically,  about 
55%,  75%  and  90%  of  total  contraction  occur  at  the  axial 
distance  of  0.2R,  0 . 5R  and  1.0R  from  the  propeller  tip  for 
all  propellers  and  for  all  loading  conditions.  Table  9 
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shows  the  contraction  ratio  with  respect  to  total  contraction 
at  the  axial  distance  of  0.2,  0.5  and  1.0  radius.  It  may  be 
also  noted  that  there  is  an  interesting  relation  between  skew 
and  ultimate  slipstream  radius,  that  is,  the  ultimate  slip- 
stream radius  is  almost  linearly  increasing  with  increasing 
skew,  as  shown  in  Figure  26.  The  relation  between  loading 
condition  and  ultimate  slipstream  radius  is  shown  in  Figure 
27. 

IV. 3.  Hydrodynamic  Pitch  and  Pitch  Angle 

With  the  LDA  system,  at  least  two  methods  are  available 
to  find  hydrodynamic  pitch  and  pitch  angle  of  tip  vortices. 
One  way  is  to  trace  the  tip  vortex  from  a selected  propeller 
blade  to  find  an  axial-angular  relation  of  that  vortex( vortex 
trace  method) . This  can  be  done  through  a series  of  plots  of 
field  point  velocity  close  to  the  slipstream  at  different 
axial  positions.  Figure  28  is  one  typical  example  describ- 
ing this  method.  The  other  method  is  to  use  the  relation 
between  axial  and  tangential  velocity  components  with  respect 
to  the  rotating  coordinate  system( velocity  method).  To  do 
this,  the  mean  velocity  components  with  respect  to  space 
fixed  coordinate  system  were  first  measured  and  converted 
to  those  with  respect  to  rotating  coordinate  system.  Both 
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methods  work  well  leading  to  the  results  with  good  agreement. 
However,  the  vortex  trace  method  is  preferred,  because  the 
result  from  the  velocity  method  tends  to  be  somewhat  dis- 
persed, while  that  of  the  vortex  trace  method  is  quite  smooth. 
Also,  the  visual  method  can  be  used  to  find  these  quantities, 
but  its  result  is  less  accurate  due  to  the  reason  explained 
in  section  IV. 2.  The  hydrodynamic  pitch  and  pitch  angle  of 
tip  vortices  obtained  from  the  vortex  trace  method  are  shown 
in  Tables  10  and  11,  and  Figures  29  through  36. 

IV. 4.  Tunnel  Free  Stream  Velocity  with  Propeller  Shaft  Housing 
and  Mean  Axial  Velocity  behind  Propeller 

Theory  assumes  no  propeller  shaft.  Actually,  there  must 
be  a shaft  and  it  generates  wake.  This  may  explain  the  reason 
why  the  axial  component  of  field  point  velocity  predicted  by 
theory  near  the  hub  vortex  is  always  considerably  higher  than 
that  of  measured,  which  will  be  shown  later.  In  order  to 
estimate  the  effect  of  wake,  tunnel  free  stream  velocities 
were  measured  with  propeller  shaft  housing  and  dummy  hub. 

They  are  shown  in  Figures  37  and  38.  Next,  the  mean  axial 
velocities  from  one  side  of  free  stream  to  the  other  side 
of  free  stream  through  the  hub  vortex  core  center  were  mea- 
sured at  several  different  axial  positions.  They  are  shown 
in  Figures  39  and  40.  Their  patterns  of  velocity  defect  are 
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comparnd  to  be  very  similar  and  it  can  be  roughly  concluded 
that  radial  region  of  0.5R  is  under  the  influence  of  the  wake. 
Also,  a spanwise  load  distribution  can  be  deduced  from  the 
mean  axial  velocity  measurement.  A comparison  between 

Figures  39  and  40  shows  the  higher  load  distribution  near  the 

* 

tip  for  highly  skewed  propellers  than  for  plain  propellers. 
This  can  be  confirmed  by  the  comparison  of  radial  circulation 
distributions  obtained  either  by  the  solutions  of  Boundary 
Value  Problem  or  by  the  tangential  velocity  measurements. 

The  methods  of  obtaining  circulation  distribution  will  be 
discussed  in  section  IV. 7. 

IV. 5.  Mean  Tangential  Velocity 

In  fact,  this  velocity  is  measured  along  the  y-axis(z=0). 
Therefore,  it  would  be  more  precise  to  state  as  "vertical 
velocity  measurement  along  the  y-axis" . However,  this  is 
equivalent  to  tangential  velocity,  particularly  to  mean 
tangential  velocity.  Figures  41  through  44  show  the  mean 
vertical  velocity  measurements  along  the  y-axis.  Their 
usefulness  will  be  shown  later.  The  mean  tangential  velocity 


* This  conclusion  applies  only  to  the  present  propeller  series, 
and  implies  that  the  present  design  method  does  not  call  for 
enough  pitch  reduction  at  the  tip. 
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and  the  instantaneous  tangential  velocities  are  almost 
identical,  that  is,  fluctuation  is  very  small  in  tangential 
velocity  except  very  close  to  tip  vortices.  Figure  45  shows 
the  comparison  of  mean  tangential  velocity  between  plain  and 
highly  skewed  propellers  measured  at  the  roll-up  point. 

IV. 6.  Roll-Up  Point 

The  position  of  roll-up  point  for  the  hub  vertex  system 
can  be  found  by  investigating  the  initial  slope  of  tangential 
velocity  through  the  hub  vortex.  Similarly,  that  of  the  tip 
vortex  system  can  be  found  by  examining  either  the  peak  axial 
velocity  or  the  fluctuation  pattern  of  axial  velocity  very 
close  to  tip  vertices.  As  shown  in  Figures  46  through  49, 
both  the  initial  slope  of  tangential  velocity  and  the  peak 
axial  velocity  at  95%  of  slipstream  radius  are  almost  line- 
arly increasing  with  increasing  axial  distance  from  the  pro- 
peller, finally  reaching  a maximum  value  at  a certain  axial 
position,  and  then  slowly  decreasing  due  to  diffusion  of 
vorticity.  The  axial  position  with  a peak  value  is  considered 
to  be  the  roll-up  point.  Thus  determined  position  of  roll-up 
points  all  fall  within  0.9R  - 1.0R  from  the  propeller  tip. 

As  mentioned  above,  another  way  for  finding  the  roll-up  point 
of  the  tip  vorticity  is  to  examine  the  fluctuation  pattern  of 
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axial  velocity  very  close  to  slipstream  radius.  As  shown 
in  Figures  50  through  53,  the  fluctuation  pattern  of  axial 
velocity  at  95%  of  slipstream  radius  first  becomes  regular 
at  a certain  axial  position,  which  is  the  same  axial  position 
as  that  of  peak  velocity.  Experience  reflects  that  the  exact 
position  of  the  roll-up  point  is  not  crucial. 


IV. 7.  Circulation 


The  tangential  velocity  measurement  is  particularly 

useful  to  find  out  the  circulation  strength  and  distribution. 

Both  initial  slope  and  definition  of  vorticity  together  with 

• * 

the  information  of  tangential  velocity  distribution  can  be 
used  to  calculate  circulation  strength.  In  (x,  r,  0)  - 
cylindrical  coordinate  system,  vorticity  is  expressed  by. 


* This  method  was  originally  suggested  by  Prof.  Widnall 
to  estimate  hub  vortex  core  size. 
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2)  there  is  no  change  in  tangential  velocity  with  change 

in  axial  distance  after  roll-up  point  i 

3)  there  is  no  change  in  radial  velocity  with  change  in 
axial  position  after  roll-up  point 

I 

< 

| 

The  first  assumption  is  true  for  mean  velocity,  and  the 
mean  velocity  was  actually  used.  The  second  and  third 

< 

assumptions  are  reasonable  for  the  portions  after  roll-up  1 

point.  In  fact,  radial  velocity  is  negligibly  small  after 

roll-up  point  except  extremely  close  to  hub  vortex  core  ! 

center.  With  the  above  assumptions,  the  vorticity  equation 

4 

is  reduced  to  the  form  of. 

= (5X,  0 ,£t)  (18) 

where , 

3X  = 4"?f(rut>  (19) 

St  = - (20) 


The  circulation  distribution  is  obtained  by  the  integration 
of  vorticity  as. 


(21) 


Thus  obtained  results  are  shown  in  Figures  54  through  57. 

In  order  to  have  some  general  aspect,  the  circulation 
strengths  have  been  calculated  by  several  different  methods 


L 


A 


7 4 


and  are  shown  in  Table  12.  Table  12  also  shows  a comparison 
of  circulation  strength  between  PUF-2  and  experiment  (vorti- 
city  integration) . Their  agreement  is  good  for  a plain  pro- 
peller, but  becomes  less  good  as  skew  increases.  This  may 
explain  the  less  good  agreement  in  propeller  thrust  between 
test  and  theory  for  highly  skewed  propellers. 

1 

IV. 8.  Field  Point  Velocity 

The  ultimate  goal  of  this  study  is  to  develop  a predic-  1 

tion  method  of  field  point  velocity  numerically  and  experi- 
mentally. In  this  regard,  the  application  of  LDA  technique 
is  quite  successful  and  the  field  point  velocity  measurement 
has  essentially  become  a routine  work.  In  fact,  most  of  pre- 
viously discussed  informations  have  been  derived  from  the 

velocity  measurements.  For  this  study,  time-averaged  field  j 

point  velocities  per  propeller  revolution  or  for  three  blade 
periods  at  over  thousand  field  points  have  been  measured  and  I 

stored  in  the  form  of  plot.  Here,  it  is  reminded  that  the 
position  of  field  points  is  expressed  following  the  ITTC  conven- 
tion shown  in  Figure  21  and  that  these  measured  velocities 
with  respect  to  space-fixed  coordinate  system  can  easily  be  con- 
verted to  those  for  propeller  fixed  (rotating)  coordinate 
system  using  the  relation  derived  in  Figure  16.  Some  examples 
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of  measurements  are  shown  in  Figures  58  through  81 . For 
convenience,  however,  most  of  measurements  around  propeller 
blade  are  presented  in  Chapter  VI  along  with  numerical  out- 
puts for  comparison.  Figures  58  through  73  show  the  typical 
fluctuation  pattern  of  axial  velocity  at  upstream  and  down 
stream  locations  rather  close  to  the  propeller  blade.  Figures 
74  through  81  show  the  particular  fluctuation  pattern  of  axial 
velocity  slightly  inside( 0 .95Rw)  and  outside( 1 . 05Rw)  of  slip- 
stream. In  fact,  these  Figures  show  the  velocity  jump  across 
the  tip  vortices.  This  velocity  jump  was  once  used  to  esti- 
mate the  tip  vortex  strength.  The  detail  method  will  not  be 
discussed  here.  It  is  simply  mentioned  that  this  method 
yields  the  value  very  close  to  that  by  initial  slope  of  mean 
tangential  velocity  method  shown  in  Table  12,  section  IV. 7. 
Also,  the  transition  from  initial  to  ultimate  fluctuation 
pattern  of  axial  velocity  slightly  inside  the  slipstream  is 
shown  in  Figures  50  through  53. 

IV. 9.  Blade  Boundary  Layer 

For  the  velocity  measurements  inside  the  slipstream, 
rather  unexpected,  but  consistent  patterns (dents)  have  been 
encountered.  Figures  82,  83  and  84  are  those  for  the  measured 
axial  velocity  clearly  showing  a dent  for  each  blade  period. 
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Although  it  seemed  more  convincing  for  them  to  be  explained 
by  viscous  effect,  there  still  remained  a possibility  that 
they  were  caused  by  vortex  systems,  because  theory  occasion- 
ally produces  dents.  Through  the  careful  investigation  of 
both  the  numerical  scheme  and  measurements,  however,  the 
conclusion  has  been  reached  that  those  are  really  velocity 
defects  caused  by  the  blade  boundary  layer.  The  dents  in 
theory  come  from  the  modeling  of  vortex  system.  In  the  pre- 
sent numerical  method,  as  will  be  discussed  in  Chapter  V, 
propeller  blades  are  replaced  by  surfaces  with  spanwise  and 
chordwise  vortex  segments,  one  with  fine  spacing  and  the 
others  with  coarse  spacing.  The  purpose  of  such  arrangements 
is  to  reduce  the  computation  time.  The  occasional  dents  in 
theory  are  due  to  this  unequal  arrangement  of  vortex  segments . 
Numerical  calculation  at  carefully  chosen  field  points  does 
not  produce  these  dents.  Furthermore,  if  the  dents  are  caused 
by  vortex  systems,  then  their  pattern  can  not  be  so  consist- 
ent as  shown  in  Figures  82,  83  and  84  with  the  variations  of 
radial  and  axial  positions.  The  dots  in  Figures  82  and  83 
are  potential  flow  approximations.  They  are  constructed 
first  by  plotting  the  theoretical  values  and  by  bringing 
down  to  match  the  measurements,  just  like  DC  voltage  change. 
This  information  of  velocity  distribution  in  the  wake  can  be 
used  to  estimate  the  blade  sectional  drag  by  applying  momentum 
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balance  for  the  properly  chosen  control  surface.  This  method 
of  determining  sectional  drag  is  not  only  very  useful,  but 
also  often  the  only  practical  way.  With  this  method,  a rough 
estimation  of  the  blade  sectional  drag  coefficient  has  been 
made  neglecting  the  contribution  from  the  static  pressure  and 
rotational  motion,  and  the  results  are  shown  in  Table  13.  As 
shown  in  Table  13,  the  estimated  sectional  drag  coefficient 
becomes  higher  toward  propeller  tip.  This  result  implies 
that  the  incorporation  of  presently  adopted  sectional  drag 
coefficient  to  the  numerical  scheme  may  not  be  fully  justi- 
fied, even  though  the  average  value  of  estimated  coefficient 
and  presently  adopted  coefficient  are  very  close  each  other. 
Here,  it  should  be  noted  that  the  estimation  is  somewhat 
affected  by  the  approximation  of  potential  flow.  Although 
those  values  are  not  expected  to  be  immediately  useful,  they 
are  at  least  valuable  to  confirm  that  the  strange  patterns 
are  really  velocity  defects  due  to  boundary  layer. 

V.  Numerical  Method 

(Numerical  Lifting  Surface  Theory) 

For  the  numerical  prediction  of  field  point  velocities. 


the  Field  Point  Velocity  Program  (FPV-9)  has  been  prepared. 
This  program  is  essentially  the  same  as  the  Propeller  Unsteady 
Force  Program  (PUF-2)  in  the  sense  of  utilizing  the  sub- 
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routines  used  in  PUF-2.  In  fact,  PUF-2  is  a consequence  of  the 
last  several  year’s  research  carried  out  by  several  people. 
Although  FPV-9  has  been  prepared  for  the  case  of  steady  flow, 
it  may  be  extended  to  unsteady  case  without  much  difficulty. 
This  program  is  composed  of  two  major  programs  - The  Boundary 
Value  Problem  Program  (BVP)  and  the  Velocity  Program.  First, 
the  propeller  geometry  is  read  in  as  an  input  in  the  Boundary 
Value  Problem  Program.  The  analytical  lifting  surface  theory 
employs  the  source  and  vortex  sheets  (continuous  distribution) 
as  a model  to  represent  propeller  blades,  and  the  trailing 
vortex  sheets  to  represent  propeller  wake.  Obviously,  such 
a modelling  is  not  practical  in  a numerical  scheme.  Therefore, 
the  continuous  distributions  of  sources  and  vortices  are  re- 
placed by  concentrated  source  and  vortex  line  segments  as 
shown  in  Figure  85.  The  vortex  spacings  on  blades  are  uniform 
both  in  the  chordwise  and  spanwise  directions.  The  roll-up 
procedures  to  the  tip  and  hub  vortices,  and  the  contraction 
of  the  slipstream  have  been  accounted  for  in  the  present  model- 
ling of  the  trailing  vortex  wake  system,  which  schematically  is 
shown  in  Figure  86.  Since  the  problem  is  basically  linear,  the 
thickness  and  lifting  problems  can  be  separated.  Applying  the 
standard  method,  the  thickness  problem  is  solved  by  distri- 
buting sources  at  each  chordwise  section.  The  source  strengths 
are  determined  by  the  mean  inflow  velocities  obtained  from  the 
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lifting  line  theory.  Thus  determined  source  strength  is 
used  exclusively  with  the  assumption  that  the  thickness 
problem  is  not  sensitive  to  the  circumferential  variation 
of  inflow  velocity.  The  lifting  problem  is  solved  utilizing 
the  existing  formulation  for  the  lifting  surface  theory, 
that  is. 

Integral  equation  for  the  unknown  vortex  strength 
Boundary  condition 
Kelvin's  condition 
Kutta  condition 

The  integral  equation  is  the  well-known  Canchy  principal  type 
singular  equation  and  particular  caution  is  required  in  the 
procedure  of  discretization.  The  boundary  condition  ( tangency 
condition)  is  satisfied  on  the  blade  camber  surface  as, 

N * V = 0 

where , 

-♦  **♦  t ->  — * 

7 = V inflow  + V source  + V vortex  (bounder  + trailer 

+ wake) 

Kelvin's  theorem  is  applied  to  locally  relate  the  trailing 
vortices  to  spanwise  vortices.  By  this  procedure,  the  number 
of  unknowns  is  reduced.  Finally,  the  Kutta  condition  is  required 
to  ensure  smooth  tangential  flow  at  the  trailing  edge.  F.om 
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the  previous  tests  and  experiences,  one-point  linear  Kutta 
condition  is  adopted.  Therefore,  the  strength  of  the  lost 
vortex  element  in  a chordwise  panel  is  linearly  dependent 
on  the  strengths  of  the  adjacent  vortices  on  the  blade  and 
in  the  wake.  This  choice  is  regarded  as  optimum,  particu- 
larly for  unsteady  cases.  The  integral  equation  is,  by 
discretization,  transformed  into  a system  of  linear  algebraic 
equations  for  the  unknown  strength  of  the  spanwise  vortices. 
Upon  solving  this  system  of  equations  and  utilizing  the  above 
conditions,  the  vorticity  distribution  over  the  blades  and 
wake  is  obtained.  In  order  to  reduce  the  computation  time 
as  briefly  mentioned  in  section  IV. 9,  the  number  of  discrete 
vortex  elements  of  the  other  blades  and  wakes  is  considerably 
reduced  as  shown  in  Figure  85 „ However,  the  calculated  velo- 
city would  not  be  noticably  affected  by  this  economization, 
particularly  for  the  field  points  close  to  a propeller  blade, 
since  the  induced  velocity  due  to  vortex  is  inversely  pro- 
portional to  the  distance  and  the  distance  between  blades  is 
sufficiently  far.  Now,  the  Velocity  Program  calculates  the 
rield  point  velocity  as  a function  of  advance  coefficient  at 
p* i f i»»d  field  points,  by  calculating  the  induced  velo- 
r'irr  ••!<•  individual  source  and  vortex  elements,  and 
-e  « j >i  • nnf  r ibutions  including  inflow  velocity. 


...  »i>  ><  .r  i . i • t * i scheme  as  in  BVP  program  is 
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applied  in  calculating  the  contributions  from  other  blades 
and  wakes.  In  summary,  the  particular  features  of  the  pre- 
sent numerical  scheme  may  be  listed  as  follows: 

1 . Incorporation  of  roll-up  procedure 

2.  Incorporation  of  contraction  of  slipstream 

3.  Boundary  condition  being  satisfied  on  the  camber 
surface  (except  thickness  effect) 

4.  One  point  linear  KuttQ.  condition 

5.  Economization 

VI . Comparison  of  Field  Point  Velocities  between 
Theory  and  Experiment 

Figures  87  through  142  show  the  field  point  velocities 
as  measured  and  as  calculated.  To  be  fair  enough,  approxi- 
mately equal  number  of  Figures  with  good  agreement  and  with 
less  good  agreement  has  been  selected  from  over  200  compar- 
isons made,  and  they  are  presented  in  order  of  upstream,  very 
close  to  leading  edge,  very  close  to  trailing  edge  and  down 
stream  for  each  four  propellers.  As  shown  in  these  Figures, 
the  agreement  at  upstream  field  is  always  excellent  for  the 
propellers  4381,  4382  and  4383  and  good  for  the  propeller 
4498.  One  exception  is  the  discrepancy  in  minimum  velocity 
as  shown  in  Figures  88,  103,  129  and  131,  as  examples.  This 
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difference  is  considered  to  come  from  both  the  imperfection 
of  theory  and  the  averaging  process  through  Educator.  Theory 
tends  to  exaggerate  the  effect  of  vortex,  particularly  when 
field  points  are  closer  to  vortex  segments,  while  slight 
change  in  propeller  RPM  caused  either  by  torsional  vibration 
of  propeller  shaft  or  by  change  of  RPM  in  driving  motor  would 
chop  off  the  sharp  change  in  velocity  through  the  averaging 

it  # , 

process  . This  phenomenon  can  actually  be  seen  if  we  visu- 
alize both  the  raw  and  averaged  signal  on  the  oscilloscope, 
as  shown  in  Figure  143.  In  Figure  143,  the  upper  figure  is 
the  raw  signal  directly  from  the  signal  processor  (tracker) 
and  the  lower  one  is  time-averaged  signal  through  the  Waveform 
Educator.  The  agreement  in  near  field  to  propeller  blade 
inside  slipstream  is  generally  good,  or  at  least  theory  does 
nicely  follow  the  experimental  pattern.  Here,  it  is  emphasized 
that  a particular  care  should  be  taken  in  making  a numerical 
prediction  of  field  point  velocities  inside  the  slipstream, 
particularly  for  the  region  from  the  trailing  edge  to  the 
roll-up  point.  As  indicated  earlier,  the  numerical  output 
is  extremely  sensitive  to  the  location  of  field  points  due 
to  the  nature  of  Cauchy  principal  type  integral  equation. 

For  the  numerical  prediction  at  arbitrary  field  points  in 


* This  phenomenon  was  previously  predicted  by  Prof.  Ezekiel. 
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this  region,  therefore,  it  is  strongly  recommended  that  a 
series  of  initial  calculations  be  made  first  at  the  geometric 
center  of  vortex  grids  and  that  proper  interpolation  be  ap- 
plied next.  To  find  out  the  position  of  geometric  center  of 
vortex  grids,  a special  program  has  been  prepared,  or  the 
Velocity  Program  will  print  out  both  Cartesian  and  polar  co- 
ordinates of  these  points  with  the  proper  indication  of 
built-in  option.  Furthermore,  experience  reflects  that  the 
position  of  initial  input  field  points  (geometric  center  of 
vortex  grids)  should  be  as  precise  as  at  least  five  decimal  places. 
This  is  one  of  limitations  of  present  numerical  scheme. 

In  fact,  most  of  velocity  measurements  inside  the  slipstream 
shown  in  the  Figures  in  this  section  were  made  at  the  imagi- 
nary geometric  center  of  vortex  grids  inside  the  tunnel  for 
direct  comparison  with  theory.  The  velocity  prediction  with 
present  numerical  scheme  is  getting  worse  with  increasing 
axial  distance  from  the  propeller  due  to  the  presently  adop- 
ted vortex  wake  model,  that  is,  the  helix  of  ultimate  wake 
is  treated  as  line  segments  instead  of  streamlines.  However, 
the  primary  interest  is  in  the  velocity  prediction  at  the 

near  field  to  propeller  blades.  F&r  field  velocity  predic- 

. . . * 
tion  is  expected  to  be  greatly  improved  when  a new  wake  model 


* This  model  is  being  prepared  by  Prof.  Kerwin 
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is  ready.  So  far,  comparisons  and  discussions  have  been 
limited  to  the  result  of  design  condition.  Even  for  off 
design  conditions,  considerable  amounts  of  measurements 
(though  much  less  than  those  for  design  condition)  were 
made  and  some  of  them  are  shown  in  Figures  144  through  157 
with  calculated  values  for  comparison.  Although  it  is 
difficult  to  make  general  comments  due  to  a limited  number 
of  comparisons,  it  seems  that  the  agreement  is  still  good 
for  the  moderately  loaded  condition( J=0 . 6) , but  it  becomes 
considerably  worse  for  the  heavily  loaded  condition( J^O . 3 ) . 
Specifically,  measurement  shows  back  flow  around  the  tip  region 
as  shown  in  Figure  147  as  an  example,  but  theory  does  not 
predict  this.  The  first  impression  is  that  present  theory 
may  not  be  adequate  for  the  prediction  of  field  point  velo- 
city at  heavily  loaded  condition. 

VII . Conclusions  and  Recommendations 

Through  this  study,  the  following  conclusions  are  made. 

1.  The  application  of  LDA  technique  is  quite  successful.  The 
previously  virtually  impossible  measurement  has  essentially 
become  a routine  work.  This  technique  can  also  be  utilized 
in  many  different  flow  measurements  as  well  as  propeller 
measurement . 
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2.  The  slipstream  radius  may  be  a function  of  many  factors 
such  as  loading  condition,  skew,  rake,  etc.  For  the  present 
series  of  propellers,  the  ultimate  slipstream  radius  almost 
linearly  increases  with  increasing  skew  or  increasing  advance 
coefficient.  With  this  study  alone,  however,  it  is  not  pos- 
sible to  derive  a generally  applicable  relationship.  This 
may  be  quite  possible  after  investigating  a wider  variety  of 
propellers  with  the  variations  including  number  of  blades, 
initial  geometry,  blade  area  ratio,  etc. 

3.  The  ultimate  hydrodynamic  pitch  and  pitch  angle  seem  to 
be  basically  determined  by  design  pitch  and  pitch  angle  and 
affected  very  little  by  skew. 

4.  For  all  loading  conditions,  most  of  contraction  of  slip- 
stream (approximately  90%)  occurs  within  the  axial  distance 
of  one  radius  from  the  propeller  plane. 

5.  Roll-up  occurs  at  the  axial  distance  of  0. 8-1.0  radius 
from  the  propeller  tip  for  design  condition.  This  axial 
position  of  roll-up  for  moderately  loaded  condition  seems 
to  be  almost  the  same  as  for  design  condition,  but  that  for 
heavily  loaded  condition  is  not  clear. 

6.  From  items  4 and  5,  it  can  be  concluded  that  there  is  a 
close  relation  between  the  position  of  roll-up  and  the  con- 
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traction  ratio  of  slipstream.  More  specifically  speaking, 
the  roll-up  points  are  approximately  the  same  as  the  positions 
of  90%  contraction  of  slipstream. 

7.  The  circulation  distribution  and  strength  obtained  by 
BVP  program  are  somewhat  affected  by  input  wake  data  such 
as  axial  position  of  roll-up,  slipstream  radius  at  roll-up 
point,  hydrodynamic  pitch  and  pitch  angle,  etc.  Tests  show 
that  these  variations  of  circulation  distribution  and 
strength  have  almost  no  effect  on  the  overall  propeller  blade 
thrust  and  torque  (output  of  PUF-2) . However,  field  point 
velocity  prediction  is  affected  as  same  degree  of  the  varia- 
tions . 

8.  The  calculated  circulation  strength  is  always  higher  than  the 
experimentally  determined  strength.  The  difference  is  small 
and  agreement  is  good  for  the  plain  propeller,  but  the  differ- 
ence becomes  greater  as  skew  increases.  This  result  may 
explain  the  less  good  agreement  in  thrust  between  theory  and 
experiment  for  highly  skewed  propellers . 

9.  The  roughly  estimated  sectional  drag  coefficient  based  on 
experimental  measurements  becomes  higher  toward  tip  instead  of 
getting  lower.  It  is  implied  by  this  result  that  the  incorpo- 
ration of  presently  adopted  sectional  drag  coefficient  to  the 
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numerical  scheme  may  not  be  fully  justified,  even  though  the 
average  value  of  estimated  coefficient  and  presently  adopted 
coefficient  used  in  the  theory  are  very  close  each  other. 

10.  For  the  field  point  velocity,  over  200  comparisons  were  , 

made  between  theory  and  experiment.  Most  of  them  are  for  the 
design  condition.  From  these  comparisons,  the  following  con- 

t 

elusion  has  been  derived:  1 

i 

For  the  design  condition,  the  agreement  is  always  excellent  for 

1 

upstream  and  outside  slipstream.  For  the  near  field  inside  the 
slipstream  including  points  very  close  to  trailing  edge,  the 
agreement  is  generally  good,  or  at  least  theory  nicely  follows 
the  experimental  pattern.  For  far  field  inside  slipstream, 
the  agreement  is  not  as  good,  but  this  region  is  not  of  much 
concern.  For  the  moderately  loaded  condition,  generally  similar 
conclusion  can  be  made  to  that  for  design  condition.  For  the 
heavily  loaded  condition,  the  agreement  is  still  not  bad  for 
upstream  and  near  field  inside  slipstream  (not  very  close  to 
trailing  edge) , but  considerably  worse  than  that  for  design 
and  moderately  loaded  conditions.  Particularly,  measurement 
shows  rather  strong  back  flow  around  tip,  but  theory  does  not 
predict  this.  The  first  impression  is  that  present  theory 
may  not  be  adequate  for  the  prediction  of  field  point  velocity 
for  heavily  loaded  condition.  However,  it  should  be  noted 
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that  the  above  conclusions  for  moderately  loaded  and  heavily 
loaded  conditions  are  derived  from  a limited  number  of  compar- 
isons and  may  be  treated  as  tentative  ones. 

11.  The  figures  of  field  point  velocity  comparison  for  near 
field  inside  slipstream  show  rather  regular  parallel  shift- 
up  of  the  numerical  results  from  measurements,  which  seems 
more  strongly  related  to  radial  position  than  to  axial  posi- 
tion. This  suggests  that  there  may  be  a general  correlation 
factor  between  measurement  and  present  theory.  With  this 
idea,  the  correlation  factor  for  the  axial  component  of  field 
point  velocity  has  been  briefly  investigated  and  shown  in 

i 

Table  14.  As  shown  in  Table  14,  the  correlation  factors  for 
two  different  propellers  at  given  radial  positions  are  remark- 
ably close( except  in  the  immediate  vicinity  of  the  hub  vortex)  « 

in  spite  of  big  difference  in  relative  axial  position  from 
propeller  blades.  The  factor  shown  in  Table  4 (as  a function 

1 

of  radial  position)  may  be  directly  useful  for  this  series  of 
propellers,  although  the  analysis  discussed  here  is  extremely 
simple . 


With  the  above  conclusions,  it  is  recommended  that  further 
research  be  conducted  in  the  following  areas « 
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1.  More  investigation  of  the  viscous  effect 

The  author  believe  that  it  is  worth  carrying  out  more 
investigation  for  the  viscous  effect.  The  viscous  effect 
seems  to  be  greater  for  lower  advance  coefficient  (higher 
loading  condition) , judging  from  blade  boundary  layer  pat- 
tern. Furthermore,  the  relation  between  blade  boundary 
layer  and  geometric  characteristics,  particularly  skew, 
should  be  carefully  studied.  Due  to  centrifugal  force, 
the  particles  tend  to  be  pushed  outward  and,  hence,  bound- 
ary layer  tends  to  be  thicker  toward  tip.  This  phenomenon 
is  particularly  important  for  highly  skewed  propellers, 
since  a considerable  part  of  the  leading  edge  may  be  already 
inside  the  boundary  layer.  The  study  of  this  effect  may  be  the 
key  to  the  solution  of  less  good  agreement  problem  in  thrust 
between  theory  and  experiment  for  highly  skewed  propellers. 

For  effective  investigation,  the  author  recommends  that  a more 
realistic  boundary  layer  pattern  be  first  constructed  by 
combining  axial  and  tangential  patterns  and  converting 
angular  scale  to  linear  scale,  and  that  comparison  of  thus 
constructed  patterns  be  made  between  plain  and  highly  skewed 
propellers . 
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2.  Correlation  study 

For  the  present  series  of  propellers,  correlation  between 
theory  and  measurement  has  been  briefly  studied  and  a correlation 
factor  has  been  derived  as  a function  of  radial  position.  It 
seems  to  be  quite  meaningful.  However,  it  is  not  expected  to 
be  generally  applicable.  The  accuracy  of  numerical  prediction 
of  field  point  velocities  will  be  improved  if  further  study  on 
this  subject  is  carried  out  to  determine  generally  applicable 
relation  between  measurement  and  present  theory. 

3.  Off-Design  Condition 

Most  of  the  measurements  and  comparisons  in  this  study  are 
made  for  the  design  condition  and  the  conclusion  is  that  the 
agreement  is  good.  However,  it  is  difficult  to  make  general 
conclusion  for  off-design  condition.  Therefore,  further 
measurements  and  comparisons  are  recommended  for  off-design 
conditions,  particularly  for  the  heavily  loaded  condition,  to 
investigate  the  applicability  of  the  present  numerical  method 
for  the  prediction  of  field  point  velocity  and  to  improve 
the  theory. 
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4.  Other  Components  of  velocity 

In  this  study,  no  comparison  has  been  made  for  other 
components  of  velocity  such  as  radial  or  tangential  compo- 
nents. The  radial  component  is  usually  considered  to  be 
small.  However,  measurement  shows  that  this  is  not  always 
true.  Even  for  the  design  condition,  almost  the  same  order  of 
free  stream  velocity  was  measured  very  close  to  blade  and  hub 
vortex.  For  the  heavily  loaded  condition,  the  radial  component 
is  usually  greater  than  free  stream  velocity  very  close  to  the 
suction  side  of  the  blades.  Although  the  agreement  between 
theory  and  measurement  for  axial  component  is  good,  that 
for  the  other  components  is  not  known.  The  comparison  of 
the  other  components  remains  to  be  done. 
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Table  1 Comparison  of  Two  Most  Common  signal  Processing 

I * 

Techniques  - Tracker  and  Counter 

p 


T 


Tracker 


Required 
S/N  ratio 

Dynamic 

Range 


Low 


2 KHZ  - 50  MHZ 


Capture 

Range 

Slew  Rate 


10%  of  full  scale 
tracking  range 


Good 


Ease  of 
Use 


Good 


Counter 


High 

(lowered  by  proper  filtering) 


1 KHZ  - 100  MHZ  if  not 
limited  by  input  filters 


not  exist  if  not  limited 
by  filters 


Limited  only  by  processing 
time  ( 1 Xis) 


Fair 


Cost 


Medium 


High 


Type 


15  mw  He-Ne  laser 


Lasef 

Wave  Length(nm) 

632, 

.8 

Initial  Beam  Diameter( D2 »mm) 

1, 

.1 

Transmitting 
( focusing ) 
lenses 

Focal  Distance(mm) 

309 

481 

Crossing  Angle(deg.) 

8.88 

5.72 

Beam  Spacing(mm) 

50 

Distance  between  Fringes(|im) 

4.087 

6.341 

Shape 

Geometric 

Diamond 

Effective 

Ellipsoid 

Beam  Diameter (mm) 

Geometric 

0.556 

0.881 

Measuring 

Effective 

0.226 

0.352 

Volume 

Length (mm) 

Geometric 

7.309 

17.652 

Effective 

2.924 

7.061 

Diameter ( mm) 

Geometric 

0.568 

0.882 

Effective 

0.227 

0.353 

Effective  Number  of  Fringes ( nearest  100) 

100 

100 

Actual  Aperture  Diameter(mm) 

0.25 

0.15 

MHZ/MPS 

0.24467 

0.15770 

MHz/fps 

0.07458 

0.04807 
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Table  4 


Tunnel  Free  Stream  Velocity  Measurement 
by  Two  Different  Methods ( ft/sec) 


Manometer 

LDA 

Manometer /LDA 

9.00 

9.049 

0.995 

11.34 

11.378 

0.997 

13.05 

13.067 

0.999 

13.73 

13.762 

0.998 

14.44 

14.376 

1.004 

16.70 

16.695 

1 .000 

17.28 

17.194 

1.005 

18.62 

18.528 

1.005 

19.66 

19.527 

1.007 

20.77 

20.647 

1.006 

100 


Table  6 NS RPC  Research  Propellers 

Principal  Design  Characteristics 
(Reproduced  from  [20"|  and  [2l")) 


Propeller  Diameter 

Hub  Diameter 

Number  of  Blades 

Expanded  Area  Ratio 

Section  Meanline 

Section  Thickness 
Distribution 


12  inches 
2.4  inches 
5 

0.725 

NACA  a=0 . 8 

NACA  66  with  NSRDC 
modified  Nose  and  Tail 


i 


r/R 

tan0^ 

C/D 

t0/c 

0.2 

1.8256 

0.174 

0.2494 

0.3 

1.3094 

0.229 

0.1562 

0.4 

1.0075 

0.275 

0.1068 

0.5 

0.8034 

0.312 

0.0768 

0.6 

0.6483 

0.337 

0.0566 

0.7 

0.5300 

0.347 

0.0421 

0.8 

0.4390 

0.334 

0.0314 

0.9 

0.3681 

0.280 

0.0239 
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4 

Table  6 (continued) 

I 


Prop. 

r/R 

0s(deg. ) 

P/D 

fQ/c 

ESI 

4.655 

1.4332 

0.0370 

ESI 

9.363 

1.4117 

0.0344 

rvj 

0.5 

13.948 

1.3613 

0.0305 

00 

co 

0.6 

18.378 

1.2854 

0.0247 

0.7 

22.747 

1.1999 

0.0199 

0.8 

27.145 

1.111. 

0.0161 

0.9 

31.575 

1.0270 

n134 

ESI 

9.29  3 

1.5124 

407 

HI 

18.816 

1.4588 

0.0385 

0.5 

27.991 

1.3860 

0.0342 

CO 

00 

co 

0.6 

36.770 

1.2958 

0.0281 

0.7 

45.453 

1.1976 

0.0230 

0.8 

54.245 

1.0959 

0.0189 

0.9 

63.102 

0.9955 

0.0159 

0.3 

9.29  3 

1.5124 

0.0407 

0.4 

18.816 

1.4588 

0.0385 

0.5 

27.991 

1.3860 

0.0342 

oo 

O' 

0.6 

36.770 

1.2958 

0.0281 

TT 

0.7 

45.453 

1.1976 

0.0230 

0.8 

54.245 

1.0959 

0.0189 

0.9 

63.102 

0.9955 

0.0159 
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Table  7 NS RPC  Research  Propellers 

I 

Some  Performance  Characteristics 
at  Design  Condition 


n.  Prop. 
\No. 

CharactV^ 

4381 

4382 

4383 

4498 

Advance 

Coeff . 

0.889 

0.889 

0.889 

0.889 

Thrust 

Loading 

Coeff. 

0.686 

0.686 

0.686 

0.686 

kt* 

0.208 

0.218 

0.231 

0.246 

iokq* 

0.419 

0.441 

0.472 

0.536 

Ip 

0.702 

0.699 

0.692 

0.649 

* NSRDC  revised  test  result 


Table  8 XSRDC  Research  Propellers 
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4498  4383  4382  4381 
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Table  9 NSRDC  Research  Propellers 

Percent  Contraction  of  Slipstream* 
at  Different  Axial  Position 


^\\^Axial  Distance 
'v  from  the 

Pro\ 

No . \ 

0.2R 

0.5R 

1 .OR 

4381 

0.322 

54.1 

74.6 

90.7 

0.600 

55.9 

75.3 

88.8 

0.889 

56.6 

75.0 

90.1 

4382 

0.285 

52.1 

75.4 

92.6 

0.600 

48.2 

72.3 

92.7 

0.889 

45.5 

71.9 

93.4 

4383 

0.252 

55.1 

77.1 

92.6 

0.600 

77.7 

93.8 

0.889 

54.6 

79.2 

94.1 

4498 

0.600 

52.8 

77.6 

93.6 

0.889 

54.5 

81.8 

96.4 

Basis  of  Percentage  : total (or  ultimate)  amount 

of  contration. 
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Table  12  NSRDC  Research  Propellers 

Comparison  of  Circulation  by  Different  Methods 


Propeller  RPM  = 900 

Advance  Coeff.(J)  = 0 . 889 ( Design) 


Prop. 

No . 

Nondimensional  Circulation(G ) x 1000 

% Dif f . * 
between 
PUF-2  and 
Exp . ( Vort . 
Integ . ) 

Theory 

Experiment 

PUF-2 

Lifting** 

Line(Lerbs 

Criterion) 

Integ . 
of 

Vorticity 

Initial 
slope  of 
Tang.  Vel . 

4381 

12.32 

10.12 

11.94 

10.44 

3.16 

4382 

12.46 

10.68 

11.98 

10.34 

4.01 

4383 

12.66 

11.24 

12.03 

11.60 

5.24 

4498 

12.77 

12.04 

12.14 

11.81 

5.19 

* Basis  of  Percentage  s Circulation  Strength  determined  by 

Experiment Vorticity  Integration 
Method ) 

**  With  measured  thrust  as  input.  Since  these  propellers 
have  unloaded  tips  relative  to  the  herbs'  optimum,  the 
latter  would  be  expected  to  result  in  a lower  value  of 
the  maximum  circulation.. 
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Table  14  NSRDC  Research  Propellers 

Comparison  of  Nondimension  Mean  Axial  Velocity 
(Correlation  between  Theory  and  Experiment) 

Advance  Coeff.(J)  = 0.889 

Speed  of  Advance(Va)  = 13.335  ft/sec 


r/R 

Prop.  4381  at  x = -0.4R 

Prop.  4383  at  x = -0.5R 

Exp. 

Theory 

Exp. 

Theory 

Exp. 

Theory 

Theory 

H 

• 

o 

0.187 

1 . 394 

0.13 

0.389 

1.399 

0.28 

0.2 

0.986 

1.392 

0.71 

1.009 

1.399 

0.72 

0.3 

1.185 

1.386 

0.86 

1.178 

1.398 

0.84 

0.4 

1.282 

1.382 

0.93 

1.258 

1.387 

0.91 

0.5 

1.328 

1.378 

0.96 

1.302 

1.385 

0.94 

0.6 

1.338 

1.374 

0.97 

1.324 

1.379 

0.96 

0.7 

1.314 

1.345 

1.324 

1.365 

0.97 

o 

• 

00 

1.252 

1.266 

1.294 

1.304 

0.99 
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Fig.  1 Basic  Concept  of  Laser-Doppler  Anemometry 


Photo 

Current 


Light  waves  scattered  from  the  illumi- 
nating beams  by  moving  particles  are  Doppler 
shifted.  The  Doppler  shift  information,  and 
hence  the  velocity  information  is  contained 
in  the  photo  current  of  the  detector  and  the 
problem  is  to  relate  the  photo  current  to 
the  particle  velocities,  beam  geometry,  etc. 


Ill 


Fig.  2 General  Relation  between  Scattered 
Light  Waves  and  Doppler  Shift 


Ref.  t Born  and  Wolf  in  1970[7l 


2 


Reference  Beam  Mode  \ 

Basic  Detection  Method  ' 


Reference  beam  method  compares  the  Doppler 
shifted  frequency  of  light  scattered  from  an  illu- 
minating beam  with  the  un  shifted  light  of  a 
reference  beam. 


fo  = fl  - f2  = l * " ?l) 

Direction  of  Measurement  i r - si 


s: 
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Fig . 4 Dual  Beam  Mode 

Dasic  Detection  Method 


Dual  beam  method  compares  the  frequencies  of  Doppler 
shifted  waves  scattered  in  the  same  direction  from  two 
different  illuminating  waves.  The  difference  or  Doppler 
frequency,  f^,  is  independent  of  the  direction  of  observation. 


fl  = fo  + X * (?  “ 

f2  = fo  + X ' (?  " 

fD  - fl  - f2  - ?•  <*2  “ 31> 

Direction  of  Measurement  i S2  - Sj  (independent  of 

the  direction 
of  observation) 

u • (sj  - s, ) = 2 ux  sin€ 
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Fig.  5 Dual  Scatter  Mode 

Basic  Detection  Method 


Dual  scatter  method  compares  the  Doppler  shifted 
frequencies  of  two  light  waves  scattered  in  different 
directions  from  the  same  illuminating  wave. 


fl  = fo  + X ’ (fl  " *0 

f2  = fo  + f ‘ " ®1) 

— ♦ 

fD  = fl  _ f2  = f ' (?1  " ?2> 
Direction  of  Measurement  i r^  - ?2 

u • ( r j - r^)  = 2 ux  sin€ 

_ 2 ux  sin€ 

fD  = X 


(Gaussian  intensity  distribution) 


Fig.  6 Characteristic0  of  Dual  Beam  Forward  Sea 
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Fig.  8 Generation  of  Fringe  Pattern 
(light  wave  interference) 


Two  interfering  » 
light  waves 


■*  _ FT-  T i[Wt  - k(xcos€  + ysin€)] 

» Ej  = -/I0  Pi  e L 

■ar  -?  i[Wt  - K(-xcosC  + ysin€)] 

^2  = VI0  P2  e 


Electric  Vector  « E = Ei  + E2 
(by  superposition) 

Intensity  1 I(x)  = (Ej  + E2)*(^l  + E2*) 

= El -El*  + E>2>E2*  + ^1*E2*  + "E2*Ei* 
= 2I0  [ 1 + "pi •P2Cos( 2kxsin€) ] 

where  K = 2rr/X 


Oscilloscope  Trace  of  Particle  Traversing 
Measuring  Volume ( Tyoical  Photodetector  Output) 


H2  cl  sinei 

Hi  C2  sin€2 

4lxl=  42X2 

X2  _ Xl/H2  _ X1 

f 2sin€2  2sin€^/42  2sin€^ 


P = total  beam  power 

h - F* 
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Fig.  12  Basic  Idea  of  Signal  Processing 
Method  of  Spectrum  Analyzer  and 
Spectrum  Broadening 


u 


Frequency 


— — « "Ideal”  signal  in  laminar  flow 

! Broadening  due  to  finite  filter  bandwidth 

» Broadening  due  to  a signal  of  finite 
length  and  noise 

Broadening  due  to  turbylence 


2 = 6^  + & * + 6^4 
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Fig.  13  Dynamic  Range  and  Capture  Bandwidth 
of  Tracker  Type  Signal  Processors 

■ 

f Velocity 


Assumption  » Ideal  system  where  system  is  always 
following  the  signal  properly (high 
slew  rate) 

A = Dynamic  Range 

(maximum  change  in  velocity  the  system  can 
follow  with  no  manual  range  change) 

B = Capture  Bandwidth 

(maximum  velocity  change  between  data  points 
that  the  system  can  measure) 
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Fig.  14  Effect  of  Noise  on  Counter 


Original  Signal 

(idealized)  


Threshold 
(after  HPF ) 


With  a given  noise  level',  accuracy  of  the 
measurement  can  be  improved  increasing 
the  number  of  cycle  time. 


125 


Fixed  Coord, 
(fixed  in  space) 


Moving  Coord. 

( f i v orl 


i rv  ktv  rv 
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Fig.  19  Thickness  Distribution  of  Tunnel 
Window  for  Laser  Experiment 
(all  dimensions  in  inches) 

4 


* 3 = 


370.209  - (1.3348)  Xx  for  Ft  = 309 

603.495  - (1.3337)  Xx  for  Ft  = 481 

with  *498( 1 .87" ) and  Correction  for  lens  mounting  case, 

Comparisons 


Ft  = 309 


Ft  = 481 


Diff . 


Meas . 

Cal 

• 

87, 

.47 

254 

• 

o 

o 

253. 

46 

107, 

.95 

225 

.82 

226. 

12 

125 

.49 

203 

.20 

202. 

71 

162, 

.08 

153 

.99 

153. 

87 

189, 

.07 

117 

00 

'T 

• 

117. 

85 

277, 

.01 

0 

.0 

0. 

45 

0.45  452.25 


*3 

Meas . 

Cal. 

456.57 

455.83 

394.59 

394.65 

332.61 

332.31 

224.41 

224.92 

144.65 

145.31 

o 

• 

o 

0.32 
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Slipstr***  Radius(Fr*ction  of  Prop.  Radius) 


Fract 


I 


of  Prop.  Radius) 
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Condition 
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Velocity  Defect  due  to  Blade  Boundary  Layer 
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Radial  Position(r)=  0.6R 
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